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I. INTRODUCTION 

Deep-inelastic scattering of leptons from nucleons has 
revealed much of what is known about quarks and glu- 
ons. The scattering of high-energy charged polarized lep- 
tons on polarized nucleons provides insight into the spin 
structure of the nucleon at the parton level. The spin- 
dependent nucleon structure functions determined from 
these measurements are fundamental properties of the 
nucleon as are the spin-independent structure functions, 
and they provide crucial information for the develop- 
ment and testing of perturbative and non-perturbative 
Quantum Chromodynamics (QCD). Examples are the 
QCD spin-dependent sum rules and calculations by lat- 
tice gauge theory. 

The first experiments on polarized electron-proton 
scattering were carried out by the E80 and E130 Col- 
laborations at SLAC W . They measured significant spin- 
dependent asymmetries in deep-inelastic electron-proton 
scattering cross sections, and their results were consis- 
tent with the EUis-Jaffe and Bjorken sum rules with 
some plausible models of proton spin structure. Sub- 
sequently, a similar experiment with a polarized muon 
beam and polarized proton target was made by the Eu- 
ropean Muon Collaboration (EMC) at CERN |||. With 
a tenfold higher beam energy as compared to that at 
SLAC, the EMC measurement covered a much larger 
kinematic range than the electron scattering experiments 
and found the violation of the EUis-Jaffe sum rule pj. 
This implies, in the framework of the Quark-Parton 
Model (QPM), that the total contribution of the quark 
spins to the proton spin is small. 

This result was a great surprise and posed a major 
problem for the QPM, particularly because of the suc- 
cess of the QPM in explaining the magnetic moments 
of hadrons in terms of three valence quarks. It stimu- 
lated a new series of polarized electron and muon nucleon 
scattering experiments which by now have achieved the 
following: 

1. inclusive scattering measurements of the spin- 
dependent structure function gf of the proton 
with improved accuracy over an enlarged kinematic 
range; 

2. evaluation of the first moment of the proton spin 
structure function, F^ — J^ g'^{x)dx, with reduced 
statistical and systematic errors; 

3. similar measurements with polarized deuteron and 
^He targets, in order to measure the neutron 
spin structure function and test the fundamental 
Bjorken sum rule for P^' — F" H; 

4. measurements of the spin-dependent structure 
function g2 for the proton and neutron; 

5. semi- inclusive measurements of final states which 
allow determination of the separate valence and sea 
quark contributions to the nucleon spin. 



The recent measurements of polarized muon-nucleon 
scattering have been done by the Spin Muon Collabo- 
ration (SMC) at CERN with polarized muon beams of 
100 GeV and 190 GcV obtained from the CERN SPS 
450 GeV proton beam and with polarized proton and 
deuteron targets. Spin-dependent cross section asymme- 
tries are measured over a wide kinematic range with rel- 
atively high Q^ and extending to low x values. The de- 
termination of gi{x, Q^) for the proton and deuteron has 
been the principal result of the SMC experiment, but (72 
and semi-inclusive measurements have also been made. 

The recent measurements of polarized electron- nucleon 
scattering have been done principally at SLAC in exper- 
iment E142 g (beam energy Ee - 19, 23, 26 GeV, ^He 
target), E143 [|,0 (beam energy Ee ~ 9, 16,29 GeV, H 
and D targets) and E154 {E^ ~ 48 GeV, ^He target). 
SLAC E155 with Ee ^ 50 GeV and polarized proton 
and deuteron targets will take data soon. The SLAC 
experiments provide inclusive measurements of gi and 
(72 over a kinematic range of relatively low Q^ and do 
not extend to very low x values. However, the electron 
scattering experiments involve very high beam intensi- 
ties and achieve excellent statistical accuracies. Hence 
the electron and the muon experiments are complemen- 
tary. Recently the HERMES experiment at DESY has 
become operational and has reported preliminary results 
with a polarized ■^He target |8|. This experiment uses a 
polarized electron beam of 27 GeV in the electron ring at 
HERA and an internal polarized gas target. Both inclu- 
sive and semi-inclusive data were obtained, and polarized 
H and D targets will be used in the future. 

In this paper, we present SMC results on the spin- 
dependent structure functions gf and (72 of the proton, 
obtained from data taken in 1993 with a polarized bu- 
tanol target. First results from these measurements were 
published in Refs. pJlCJ]. We use here the same data sam- 
ple but present a more refined analysis; in particular, we 
allow for a Q^-evolution of the gi structure function as 
predicted by perturbative QCD. SMC has also published 



results on the deuteron structure function g^ [|11[-|13[ and 
on a measurement of semi-inclusive cross section asym- 
metries p4i. For a test of the Bjorken sum rule we refer 
to our measurement of gf . 

The paper is organized as follows. In Section O, we re- 
view the theoretical background. The experimental set- 
up and the data-taking procedure are described in Sec- 
tion 



HI 



In Section IV we discuss the analysis of cross 
section asymmetries and in Section |V| we give the eval- 
uation of the spin-dependent structure function g^ and 
its first moment. T he r esults for ^2 ^re discussed in Sec- 
tion VI . In Section VII we combine proton and deuteron 
results to determine the structure function gf of the neu- 
tron and to test the Bjorken sum rule. In Section VIII we 



interpret our results in terms of the spin structure of the 
proton. Finally, we present our conclusions in Section IX . 



II. THEORETICAL OVERVIEW 

A. The cross sections for polarized lepton-nucleon 
scattering 

The polarized deep-inelastic lepton-nucleon inclusive 
scattering cross section in the one-photon exchange ap- 
proximation can be written as the sum of a spin- 
independent term a and a spin-dependent term Act and 
involves the lepton helicity hi = ±1: 



a — a hf/S.a. 

2 



(2.1) 



For longitudinally polarized leptons the spin S^ is along 
the lepton momentum k. The spin-independent cross 
section for parity-conserving interactions can be ex- 
pressed in terms of two unpolarized structure functions 
Fi and ¥2- These functions depend on the four mo- 
mentum transfer squared Q^ and the scaling variable 
X = Q^ /2Mv, where v is the energy of the exchanged 
virtual photon, and A/ is the nucleon mass. The double 
differential cross section can be written as a function of 
X and Q2 PI: 



d^a 



Ana^ 



dxdQ^ Q^x 



xy'{l-^)F,ix,Q^) + 



{l-y-^)F2{x,Q') 



(2.2) 



where mg is the lepton mass, 
system, and 

2Mx 



y = v / E in the laboratory 



(2.3) 



The spin-dependent part of the cross section can be 
written in terms of two structure functions 51 and (72 
which describe the interaction of lepton and hadron cur- 
rents. When the lepton spin and the nucleon spin form 
an angle -0, it can be expressed as pi] 



Aa = cos ij} Acr|| -I- sin ^ cos (j) ^cr± , 



(2.4) 



where is the azimuthal angle between the scattering 
plane and the spin plane (Fig. m). 

The cross sections Acr|| and Acr^ refer to the two con- 
figurations where the nucleon spin is (anti)parallel or or- 
thogonal to the lepton spin; Aa\\ is the difference be- 
tween the cross sections for antiparallel and parallel spin 
orientations and Acrj^ — — /i^ Actt/cos^, the difference 
between the cross sections at angles and (/) + tt. The 
corresponding differential cross sections are given by 
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Spin Plane 




Scattering Plane 

FIG. 1. Lepton and nucleon kinematic variables in polar- 
ized lepton scattering on a fixed polarized nucleon target. 
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(2.6) 



For a high beam energy E, j is small since either x is 
small or Q^ large. The structure function gi is therefore 
best measured in the (anti)parallel configuration where 
it dominates the spin-dependent cross section; (72 is best 
obtained from a measurement in the orthogonal config- 
uration, combined with a measurement of gi. In all 
formulae used in this article we consider only the sin- 
gle virtual-photon exchange. The interference effects be- 
tween virtual Z'^ and photon exchange in deep-inelastic 
muon scattering have been measured [|l7| and found to 
be small and compatible with the standard model expec- 
tations. They can be neglected in the kinematic range of 
current experiments. 



B. The cross section asymmetries 

Th e spin-dependent cross section terms, Eqs. ( ^.Sj ) and 
( ^.6|) , make only a small contribution to the total deep- 
inelastic scattering cross section and furthermore their 
contribution is, in general, reduced by incomplete beam 
and target polarizations. Therefore they can best be de- 
termined from measurements of cross section asymme- 
tries in which the spin-independent contribution cancels. 
The relevant asymmetries are 



^11 



Ao^ 
2a 



A, 



Aa± 
~2d 



(2.7) 



which are related to the virtual photon-proton asymme- 
tries A\ and A^ by 



where 



A^ = d{A2~iA^), (2., 



A2 = 



'^1/2 - cr3/2 _ ,91 - 7 52 



0'l/2 + 0-3/2 
0"l/2 + 03/2 



= 7- 



.9 1 +ff2 



(2.9) 



In Eqs. ( |2.S| ) and (p.9|), D is the depolarization factor of 
the virtual photon defined below and d, 77 and f are the 
kinematic factors: 



1 - y/2 

7(l-y-72y2/4) 



7? 



(l-2//2)(l + 72y/2)^ 

7(1 -y/2) 



I+7V2 



(2.10) 
(2.11) 
(2.12) 



The cross sections cri/2 and 173/2 refer to the absorption 
of a transversely polarized virtual photon by a polarized 
proton for total photon-proton angular momentum com- 
ponent along the virtual photon axis of 1/2 and 3/2, re- 
spectively; a^^ is an interference cross section due to the 
helicity spin-flip amplitude in forward Compton scatter- 
ing [|l8[. The depolarization factor D depends on y and 
on the ratio R = ai^/aT of longitudinal and transverse 
photoabsorption cross sections: 



£» = 



y(2-t/)(l+7V2) 



2/2(1 +y)(l-2m2/Q2)+2(l-y- 



-7V/4)(l + i?) ■ 
(2.13) 



From Eqs. ( p.8[ ) and (^), we can express the virtual 
photon-proton asymmetry Ai in terms of gi and A2 and 
find the following relation for the longitudinal asymme- 
try: 



A 



7^ = " 






(2.14) 



The virtual-photon asymmetries are bounded by positiv- 
ity relations \Ai\ < 1 and 1^2 1 < /R [|l||. When 



the term proportional to A2 is neglected in Eqs. (2.?) 
and ( ^.14 ), the longitudinal asymmetry is related to Ai 
and gi by 



Ai 



D 
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An 



Fi 1 -I- 72 D ' 



(2.15) 



respectively, where Fi is usually expressed in terms of F2 
and R: 



_ 1 + 7^ 
2x{l + R) 



(2.16) 



These relations are used in the present analysis for the 
evaluation of gi in bins of x and Q^, starting from the 
asymmetries measured in the parallel spin configuration 
and using parametrizations of F2{x, Q^) and R{x, Q^)- 



The virtual photon-proton asymmetry A2 is evaluated 
from the measured transverse and longitudinal asymme- 
tries A\\ and A±: 
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r/e V d 



.A|| 

■u 



(2.17) 



From Eqs. ( |2.3[ ) and (p.9[), A2 has an explicit 1/vQ^ 
dependence and is therefore expected to be small at high 
energies. The structure function Q2 is o btain ed fro m the 
measured asymmetries using Eqs. ( ^.9|) and ( 2.17 ). 



C. The spin-dependent structure function gi 

The significance of the spin-dependent structure func- 
tion gi can be understood from t he v irtual photon asym- 
metry Ai. As shown in Eq. ( |2.9| ), Ai ~ gi/Fi, or 
fi/2 — cr3/2 oc gi. In order to conserve angular momen- 
tum, a virtual photon with helicity -1-1 or —1 can only 
be absorbed by a quark with a spin projection of — ^ 
or -l-i, respectively, if the quarks have no orbital angu- 
lar momentum. Hence, gi contains information on the 
quark spin orientations with respect to the proton spin 
direction. 

In the simplest Quark-Parton Model (QPM) , the quark 
densities depend only on the momentum fraction x car- 
ried by the quark, and gi is given by 



gi{x) 



IT."' 



CiAq.ix), 



(2.18) 



where 



Aq,{x) = q+ix) - qi{x) + q+(x) - qf{x), (2.19) 

qt{qi^) and q^{qi~) are the distribution functions of 
quarks (antiquarks) with spin parallel and antiparallel 
to the nucleon spin, respectively, e^ is the electric charge 
of the quarks of flavor i; and Uf is the number of quark 
flavors involved. 

In QCD, quarks interact by gluon exchange which gives 
rise to a weak Q^ dependence of the structure functions. 
The treatment of gi in perturbative QCD follows closely 
that of unpolarized parton distributions and structure 
functions Q- At a given scale Q"^, gi is related to the po- 
larized quark and gluon distributions by coefficient func- 
tions Cq and Cg through Q 



"-f p2 



'dy 



g,ix,t)^^Y.- / - CS(|,a,(i))AI](y,i) 

fc=i ^f -'^ y 



+ 2nfCgi^,a,it))Ag{y,t) 
+ C^^i^,a,it))Aq^^iy,t) 



(2.20) 



In this equation, t = ln(Q2/A2), as is the strong coupling 
constant, and A is the scale parameter of QCD. The su- 
perscripts S and NS, respectively, indicate flavor-singlet 



and non-singlet parton distributions and coefRcient func- 
tions; A(7(x, t) is the polarized gluon distribution and AE 
and Aq^^ are the singlet and non-singlet combinations 
of the polarized quark and antiquark distributions 



(2.21) 




Aq'^^{x,t) 



Aq,{x,t). (2.22) 



The t dependence of the polarized quark and gluon 
distributions follows the Gribov-Lipatov-Altarelli-Parisi 
(GLAP) equations |2^p3]. As for the unpolarized distri- 
butions, the polarized singlet and gluon distributions are 
coupled by 



Pgg in 1^,^. This formalism allows a complete Next- 
to-Leading Order (NLO) QCD analysis of the scaling vi- 
olations of spin-dependent structure functions. 

In QCD, the ratio gi/Fi is Q^-dependent because the 
splitting functions, with the exception of Pqq, are differ- 
ent for polarized and unpolarized parton distributions. 
Both Pgq and Pgg are different in the two cases be- 
cause of a soft gluon singularity at x = which is only 
present in the unpolarized case. However, in kinematic 
regions dominated by valence quarks, the Q^ dependence 
of gi/Fi is expected to be small [p6[. 



D. The small-x behavior of gi 

The most important theoretical predictions for polar- 
ized deep-inelastic scattering are the srmi rules for the 
nucleon structure functions gi. The evaluation of the 
first moment of gi , 



di "• ' ' ^TT J, y 



P|,(f,a,(i))AS(2;,i) 



qq^y 



-Ag{x,t) = 

at ZTT 



+ 2nfPggi^,as{t))Ag{y,t) 



(2.23) 



y 



Pgg{^,aAt))AJ:iy,t) 
+ Pgg{^,a,{t))Ag{y,t) 



(2.24) 



whereas the non-singlet distribution evolves indepen- 
dently of the singlet and gluon distributions: 



±Aq^^ix,t) ^^ J^ ^P^^^i^-,a,it))Aq^^{y,t). 

(2.25) 

Here, Pij are the QCD splitting functions for polarized 
parton distributions. 



Expressions ( |2.20D , ( p. 23]) , ( |2.24| ) and ( |2.25| ) are valid in 
all orders of perturbative QCD. The quark and gluon dis- 
tributions, coefRcient functions, and splitting functions 
depend on the mass factorization scale and on the renor- 
malization scale; we adopt here the simplest choice, set- 
ting both scales equal to Q^. At leading order, the coef- 
ficient functions are 



Cf(|,a,) = CO-^S(^^„^) 



^(1-f) 



g y y ^ -^ 







(2.26) 



Note that gi decouples from Ag in this scheme. 

Beyond leading order, the coefficient functions and the 
splitting functions are not uniquely defined; they depend 
on the renormalization scheme. The complete set of co- 
efficient functions has been computed in the MS renor- 
malization scheme up to order a^ |2^. The 0{a^) cor- 
rections to the polarized splittini 
have been computed in Ref. 123 



functions Pgg and Pgg 
and those to Pgg and 



ri(Q') 



gi{x,Q'^)dx, 



(2.27) 



requires knowledge of gi over the entire x region. Since 
the experimentally accessible x range is limited, extrap- 
olations to X = and x = 1 are unavoidable. The lat- 
ter is not critical because it is constrained by the bound 
1^1 1 < 1 and gives only a small contribution to the inte- 
gral. However, the small- a; behavior of gi{x) is theoret- 
ically not well established and evaluation of Fi depends 
critically on the assumption made for this extrapolation. 

From the Regge model it is expected that for Q^ <^ 
2Mi/, i.e. x ^ 0, gf+g'^ and g^-g'^ behave like x-"" [^, 
where a is the intercept of the lowest contributing Regge 
trajectories. These trajectories are those of the pseu- 
dovector mesons /i for the isosinglet combination, g^+gi 
and of fli for the isotriplet combination, g^ — g", respec- 
tively. Their intercepts are negative and assumed to be 
equal, and in the range —0.5 < a < 0. Such behavior has 
been assumed in most analyses. 

A flavor singlet contribution to gi{x) that varies as 
(2 In J — 1) |28) was obtained from a model where an ex- 
change of two nonperturbative gluons is assumed. Even 
very divergent dependences like gi{x) ex (a; In a;)""'^ were 
considered |£9[. Such dependences are not necessarily 
consistent with the QCD evolution equations. 

Expectations based on QCD calculations for the be- 
havior at small- a; of gi{x, Q"^) are two-fold: 

• resummation of standard Altarelli-Parisi correc- 



tions gives 30-32 



gi{x,Q'^) ~exp 



A^\n{a,{Ql)la,{Q^))Ml/x) 

(2.28) 
for the non-singlet and singlet parts of gi . 



resummation of leading powers of In l/x gives 






(2.29) 
(2.30) 



for the non-singlet ||33| and singlet |^ parts, re- 
spectively. 



E. Sum-rule predictions 

1. The first moment of gi and the Operator Product 
Expansion 

A powerful tool to study moments of structure func- 
tions is provided by the Operator Product Expan- 
sion (OPE), where the product of the leptonic and 
the hadronic tensors describing polarized deep-inelastic 
lepton-nucleon scattering reduces to the expansion of 
the product of two electromagnetic currents. At lead- 
ing twist, the only gauge- invariant contributions are due 
to the non-singlet and singlet axial currents |g^,Q. If 
only the contributions from the three lightest quark fla- 
vors are considered, the axial current operator Ak can be 
expressed in terms of the SU(3) flavor matrices A^ (fc = 
1,...,8) and Ao = 2/ as pi 



and the first moment of gi is given by 



(2.31) 



vr?'"'(g 
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c^HQ') r 
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+ (-) {ps\ Al \ps) 



^ {ps\ Al \ps) 



(2.32) 



where Cf^ and Cf are the non-singlet and singlet coeffi- 
cient functions, respectively. The proton matrix elements 
for momentum p and spin s, {ps\ A^ |ps), can be related 
to those of the neutron by assuming isospin symmetry. In 
terms of the axial charge matrix element (axial coupling) 
for flavor qi and the covariant spin vector s^, 



St,a,{Q'^) = (j>s\ q^'y5'l^J.q^ \ps) 
they can be written as 



{ps\Al\ps) =^a3 = ^(a„-ad) = ^ 

{ps\ Al \ps) ^ ^«8 = ^("" + ad - 

{ps\ Al \ps) = Sf^aa = Sp(au + a^ + as), 
= s^ao(Q^), 



gA 



9v 

2as), 



(2.33) 

(2.34) 

(2.35) 
(2.36) 



where the Q^ dependence of a„, ad a nd a s is implied 
from now on and is discussed in Section II F . The matrix 
element a^ in Eq. ( ^.34 ) under isospin symmetry is equal 
to the neutron /?-decay constant qa/ 9v- If exact SU(3) 
symmetry is assumed for the axial-flav or oc tet cu rrent , 
the axial couplings 03 and ag in Eqs. ( 2.34 ) and ( 2.35| ) 
can be expressed in terms of coupling constants F and 
D, obtained from neutron and hyperon /3-decays H, as 



03 



F + D as^3F-D. 



(2.37) 



The effects of a possibl e SU(3 ) symmetry breaking will 
be discussed in Section [VIII B| . 

The first moment of the polarized quark distribution 
for flavor qi, that is Ag^ = J Aqi{x)dx, is the contribution 
of flavor qi to the spin of the nucleon. In the QPM ai is 
interpreted as Ag^ and oq as AE = Au -I- Ad -I- As. In 
this framework, the moments Ou, fld, Os • ■ • are bound by 
a positivity limit given by the corresponding moments of 
M, d, s, . . . obtained from unpolarized structure functions. 
In Section [I F we will see that the U(l) anomaly modifies 
this simple interpretation of the axial couplings. 

When Q^ is above the charm threshold {2mc)^, four 
flavors must be considered and an additional proton ma- 
trix element must be defined: 



{ps\ Al^ \ps) = ^-/|(au + ad + as~ 3ac) = — ^ 



ai5. 



(2.38) 



while the singlet matrix element becomes Sf^{a^ + Od 
Os -I- ac). 



2. The Bjorken sum rule 

The B jorke n sum rule B is an immediate consequence 
of Eqs. (IH) and ( |l3^ ). In the QPM where Cf^ = 1, 



pP pn 
^1 ^1 



6 



9v 



(2.39) 



In this form, the sum rule was first derived by Bjorken 
from current algebra and isospin symmetry, and has since 
been recognized as a cornerstone of the QPM. 

The Bjorken sum rule is a rigorous prediction of QCD 
in the limit of infinite momentum transfer. It is subject to 
QCD radiative corrections at finite values of Q^ |35| , P7| . 
These QCD corrections have recently been computed up 
to ©(a?) [Q and the 0{a^) correction has been esti- 
mated [ p9[ . Since the Bjorken sum rule is a pure flavor 
non-singlet expression, these corrections are given by the 
non-singlet coefficient function C^^: 



T\ - Tl = 



1 



9A 



9v 



cl 



NS 



(2.40) 



Beyond leading order, C^^ depends on the number of 
fiavors and on the renormalization scheme. Table | shows 
the coefficients cf ^ of the expansion 



TABLE I. Higher-order coefficients of the non-singlet and singlet coefficient functions Cf^ and Cf in the MS scheme. The 
coefficie nts c^ and cf are estimates; cf is unknown for n/ = 4 flavors. The quantities ag" and ao{Q^) are discussed in 
Section III E A. 
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cfs 


non- 


singlet 


cr 


cf 


singlet (a^) 
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cf 


cf 


singlet (ao(Q^)) 
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cf 


3 

4 


1.0 
1.0 


3.5833 
3.2500 


20.2153 
13.8503 


130 

68 


0.3333 
0.0400 


0.5496 
-1.0815 
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1 
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1.0959 
-0.0666 
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in the MS scheme. 
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NS\ 



as{Q^ 



(2.41) 



Since 00 = 03 + 3os, the assumption Og = is equivalent 
to ag = as = 3F—D. The quantity 3F—D is independent 
of Q^, so the assumption oq = og should be made for 
fljf ^ ao(Q^ = 00) @ 0. The coefficients cf in the 
third column of TableTshould be used to comp ute the 
coefficient function Cf that appears in Eq. ( ^.44| ). 



3. The Ellis-Jaffe sum rules 

In the QPM the coefficient functions are equa l to unity 
and assu ming exact SU(3) symmetry (Eq. ( 2.37 )) the ex- 
pression (2.32) can be written: 



^p(n) _ 



-+H^(^- 



D)^y,(^F 



D) + -o, . (2.42) 



This relation was derived by Ellis and Jaffe |g] . With the 
additional assumption that Og = 0, which in the QPM 
means As = 0, they obtained numerical predictions for 
rP and r'j". The EMC measurement ||] showed that F^* is 
smaller than their prediction which in the QPM implied 
that AS, the contribution of quark spins to the proton 
spin, is small. This result is at the origin of the current 
interest in polarized deep-inelastic scattering. 

The moments of gi and the Ellis-Jaffe predictions are 
also subject to QCD radiative corrections. The coef- 



ficient function Cf^ (Eq. (2.41)) used for the Bjorken 
sum rule also applies to the non-singlet part. The addi- 
tional c oeffic ient function Cf for the singlet contribution 
in Eq. (2.32) has been computed up to O^aj.) |3^ and 
the Ojaj) term has also been estimated for nj — i fla- 
vors lEoi: 



4- Higher twist effects 

As for unpolarized structure functions, spin-dependent 
structure functions measured at small Q^ are subject to 
higher twist (HT) effects due to nonperturbative contri- 
butions to the lepton-nucleon cross section. In the anal- 
ysis of moments and for not too low Q^, such effects are 
expressed as a power series in l/Q^: 



rr = ia(") + i^(o(^)+4d(^)+4/(^)) + o(|J 



1 



,(0) 



HT. 



(2.45) 



Here 0^"'^^ S-"^^ and Z^^-* are the reduced matrix ele- 
ments of the twist-2, twist-3 and twist-4 components, 
respectively, and M is the nucleon mass. The values of 
o^^) and S'^' for proton and deuteron have recently been 
measured [BI| from the second moment of 51 and (72, and 
found to be consistent with zero. Several authors have es- 
timated the HT effects for Fi [||-||| and for the Bjorken 
sum rule |l45| , [46| . In the literature, there is a consensus 
that such effects are probably negligible in the kinematic 
range of the data used to evaluate Fi in this paper. 



cl 
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(2.43) 



and the coefficients cf are shown in Table pi The QCD- 
corrected Ellis-Jaffe predictions for Ug = become 



-^p(n) 



Cf 



NS 



+ (-) 
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36 



(3F - D) 



--Cf{3F-D). (2.44) 



F. The physical interpretation of ao and the U(l) 
anomaly 

In the simplest approximation, the axial coupling 
o-oiQ^) is expected to be equal to AE, the contribu- 
tion of the quark spin to the nucleon spin. However, 



*In Ref. tefl, a^ and ao{Q ) are referred to as Ei„i, 
E((5^), respectively. 



and 



in QCD the U(l) anomaly causes a gluon contribution 
to ao{Q^) Wn¥^ as well which makes AS dependent on 
the factorization scheme, while ag is not. The total frac- 
tion of the nucleon spin carried by quarks is the sum of 
AS and Lq, where Lq is the contribution of quark orbital 
angular momentum to the nucleon spin. Recently, it was 
pointed out [ pO[ that this sum is scheme-independent be- 
cause of an exact compensation between the anomalous 
contribution to AS and to Lq. 

The decomposition of gq into AS and a gluon contri- 
bution is scheme-dependent ||5l|. In the Adler-Bardeen 
(AB) ||5|] factorization scheme~||53[| 



ao(Q^) = AS - n/ 



2tt 



Ag{Q' 



(2.46) 



where the last term was originally identified as the 
anomalous gluon contribution |47|-|49t|. In this scheme 
AS is independent of Q^; however it cannot be obtained 
from the measured oq without an input value for Ag. In 
other schemes AS is equal to ao((3^) but then it depends 
on Q^ [|l|. The differences between these two schemes 
do not vanish when Q^ — > oo because as{Q^)Ag remains 
finite when Q^ ^ oo @ . 



G. The spin-dependent structure function (72 

Phenomenologically, the structure function 52 can be 
understood from the spin-flip amplitude that gives rise 
to the interference asymmetry A2 (x gi + 52 of Eq. (2.9), 
owing to the absorption of a longitudinally polarized pho- 
ton by the nucleon. There are two mechanisms by which 
this can occur [ p4| . In the first, allowed in perturba- 
tive QCD, the photon is absorbed by a quark, causing 
its helicity to flip, but since helicity is conserved for 
massless fermions, this process is strongly suppressed for 
small quark masses. In the second, which is of a non- 
perturbative nature, the photon is absorbed by coherent 
parton scattering where the final-state quark conserves 
helicity by absorption of a helicity —1 gluon. 

Wandzura and Wilczek have shown |55|| that 52 can be 
decomposed as 



g2{x,Q^)^gY'^{x,Q^) + -92{^,Q^). 
The term g^^ is a linear function of gi, 



(2.47) 



g7^{x,Q^)^-gi{x,Q^) + j gi{t.Q^)\. (2.48) 

The term g2 is due to a twist-3 contribution in the 
OPE ||l^ and is a measure of quark-gluon correlations 
in the nucleon ||56| . 

In the simplest QPM, 32 vanishes because the masses 
and transverse momenta of quarks are neglected. The 
predictions of improved quark-parton models which take 
these aspects into account depend critically on the as- 
sumptions made for the quark masses and the nucleon 
wave function IM . 



The Burkhardt-Cottingham sum rule predicts that the 
first moment of (72 vanishes for both the proton and the 
neutron IBTI: 



r. 



g2{x)dx = 



(2.49) 



This sum rule is derived in Regge theory and relies on 
assumptions that are not well established. Its validity 
has therefore been the subject of much debate in the 
recent theoretical literature |16 58,M. 



III. EXPERIMENTAL METHOD 
A. Overview 

The experiment involves principally the measurement 
of cross section asymmetries for inclusive scattering of 
longitudinally polarized muons from polarized protons 
in a solid butanol target (Fig. gI). The energy of the 
incoming positive muons, 190 GeV, is measured with 
a magnetic spectrometer in the Beam Momentum Sta- 
tion (BMS). The scattered muons are detected in the 
Forward Spectrometer (FS). They are identified by coin- 
cident hits in arrays of hodoscopes located upstream and 
downstream of a hadron absorber; their momenta are 
measured with a large-acceptance, high-resolution mag- 
netic spectrometer. The beam polarization is measured 
with a polarimeter located downstream of the FS. The 
high energy of the beam provides a kinematic coverage 
down to x ~ 0.003 for Q^ > 1 GeV^, and a high aver- 
age Q^ . A small data sample was collected with a beam 
energy of 100 GeV and transverse target polarization for 
the measurement of the asymmetry A^ . 

The counting-rate asymmetries measured in this ex- 
periment vary from 0.001 to 0.05 depending on the kine- 
matic region. To assure that the asymmetries measured 
do not depend on the incident muon flux, the polarized 
target is subdivided into two cells which are polarized 
in opposite directions. Frequent reversals of the target 
spin directions in both cells strongly reduce systematic 
errors arising from time-dependent variations of the de- 
tector efficiencies. Such errors are further reduced by the 
high redundancy of detectors in the forward spectrome- 
ter. The muon beam polarization is not reversed in this 
experiment. 

The statistical errors of the counting-rate asymmetries 
are proportional to {Pf^Pt)^^ {N)~^^^ , where P^ and Pt 
are the beam and target polarizations, respectively, and 
N is the number of events. Hence high values of P^ and 
Pt as well as high N are important. 



B. The muon beam 

The SMC experiment (NA47) is installed in the up- 
graded muon beam M2 of the CERN SPS ii. A beryl- 



lium target is bombarded with 450 GeV protons from 
the SPS and secondary pions and kaons are momentum- 
selected and transported through a 600 m long decay 
channel where for 200 GeV about 5 percent decay 
into muons and neutrinos. The remaining hadrons 
are stopped in a 9.9 m long beryllium absorber for the 
190 GeV muon beam. Downstream of the absorber, 
muons are momentum selected and transported into the 
experimental hall. 

The beam intensity was 4 x 10^ muons per SPS pulse; 
these pulses are 2.4 s long with a repetition period of 
14.4 s. The beam spot on the target was approximately 
circular with a r.m.s. radius of 1.6 cm and a r.m.s. mo- 
mentum width of « 2.5%. The momentum of the in- 
cident muons is measured for each trigger in the BMS 
located upstream of the experimental hall (Fig. 0). The 
BMS employs a set of quadrupoles (Q) and a dipole (B6) 
in the beam line, with a nominal vertical deflection of 
33.7 mrad. Four planes of fast scintillator arrays (HB) 
upstream and downstream of this magnet are used to 
measure the muon tracks. The resolution of the momen- 
tum measurement is better than 0.5%. 

The beam is naturally polarized because of parity vi- 
olation in the weak decays of the parent hadrons. For 
monochromatic muon and hadron beams, the polariza- 
tion is a function of the ratio of muon and hadron ener- 
gies i|: 



P,=±- 



'tt.K 



(1 



2B„ 



-)ml 



(3.1) 



'tt.K 



where the — and -I- signs refer to positive and negative 
muons, respectively (Fig. |^). For a given pion energy, 
the muon intensity depends on the ratio Et^^k/E^; this 
ratio was optimized using Monte Carlo simulations of the 
beam transport [|62| , |63| to obtain the best combination of 
beam polarization and intensity. 



C. Measurement of the beam polarization 

A polarimeter downstream of the muon spectrometer 
allows us to determine the beam polarization by two dif- 
ferent methods. The first involves measuring the en- 
ergy spectrum of positrons from muon decay in flight, 
/i"*" — > e^DnVc, which depends on the parent-muon polar- 
ization [|64| . The second method involves measuring the 
spin-dependent cross section asymmetry for elastic scat- 
tering of polarized muons on polarized electrons Bq] . The 
two methods require different layouts for the polarimeter, 
and thus cannot be run simultaneously. 



1. Polarized-muon decay 

The energy spectrum of positrons from the decay /i"*" — > 
b^VqV^ pql can be expressed in terms of the ratio of 
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FIG. 2. Schematic layout of the muon beam and forward 
spectrometer. The individual detectors are discussed in the 
text (see Table ill). In (b), Bll is a compensating dipole that 



is used only when taking data with transverse target polar- 
ization. In (c), B8 is the forward spectrometer magnet and 
referred to as the FSM in the text. A right-handed coordi- 
nate system is used with its origin at the center of B8. The 
X-axis points along the beam direction, and the z-axis points 
upwards (out of the page in (b) and (c)). 
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FIG. 3. Muon polarization P^j as a function of muon beam 
energy E^ Ml for a monochromatic pion beam of 205 GeV 
(solid line)(Eq. (3.1)), and mean P^ vs. E^ as calculated by 
beam transport simulations [M (dashed line). 



polarimeter spectrometer matches the energy deposition 
in the LGC. 

The measured positron spectrum is corrected for the 
overall detector response. The response function is ob- 
tained from a Monte Carlo simulation that generates 
muons according to the measured beam phase space. The 
simulation accounts for radiative effects at the vertex and 
external bremsstrahlung, the geometry of the set-up, and 
chamber efficiencies. The Monte Carlo events were pro- 
cessed using the same procedure applied to the real data. 
The response function is obtained by dividing th e M onte 
Carlo spectrum by the Michel spectrum of Eq. (3.2). 



The polarization P^ can be determined by fitting 



Eq. ( |3.2| ) to the measured decay spectrum corrected for 
the detector response. Figure ^ shows the sensitivity of 
the Michel spectrum to the muon polarization. The sys- 
tematic error in the P^ determination is mainly due to 
uncertainties in the response function, the main contribu- 
tions to which arc uncertainties in the MWPC efficiencies 
and in the background rejection. Background due to ex- 
ternal 7-conversion, /i+ — > /i^7 — > /i"^e"^e~, is measured 



where A^o is the number of muon decays. 

The polarimeter configuration for this measurement is 
shown in Fig. H(a). It consists of a 30 m long evacu- 
ated decay volume, followed by a magnetic spectrome- 
ter and an electromagnetic calorimeter to measure and 
identify the decay positrons. The beginning of the de- 
cay path is defined by the shower veto detector (SVD) 
which consists of a lead foil followed by two scintillator 
hodoscopes. Along the decay path, tracks are measured 
with multiwire proportional chambers (MWPC). The de- 
cay positrons are momentum analyzed in a magnetic 
spectrometer consisting of a 6 meter long small-aperture 
dipole magnet followed by another set of MWPC. This 
spectrometer and the BMS, which measures the parent 
muon momentum, were intercalibrated in dedicated runs 
to 0.2%. A lead glass calorimeter (LGC) is used to iden- 
tify the decay positrons. 

The trigger requires a hit in each SVD plane, in co- 
incidence with a signal from the LGC above a threshold 
of about 15 GeV. Events with two or more hits in both 
planes within a 50 ns time window are rejected. This 
suppresses background from incident positrons originat- 
ing upstream of the polarimeter and rejects events with 
more than one muon. 

In the off-line analysis, events whose energy E^ was 
measured in the BMS and experienced a large energy 
loss in the SVD are rejected. A single track is required, 
both upstream and downstream of the magnet. To reject 
muon decays inside the magnetic field volume, the up- 
stream and downstream tracks are required to intersect 
in the center of the magnet. Decay positrons are iden- 
tified by requiring that the momentum measured by the 
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FIG. 4. Schematic layout of the beam polarimeter for 
the muon decay measurement (a) and for the muon-electron 
scattering measurement (b). The different components of the 
apparatus are discussed in the text. The lead glass electro- 
magnetic calorimeter and the shower veto detector are labeled 
as LGC and SVD, respectively. 
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using the charge-conjugate process with a fi" beam and 
was found to be negHgible. Other contributions to the 
systematic error arise from uncertainties in ?/e, in radia- 
tive effects at the vertex and in the ahgnment of the wire 
chambers. 



2. Polarized muon-electron scattering 

In QED at first order, the differential cross section for 
elastic scattering of longitudinally polarized muons off 
longitudinally polarized electrons is [p9| 



da 
dyf,e 
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(3.3) 
the classical electron 



(1 + ml/2m,E^)-^ 



where 171^ is the electron mass, r, 
radius, j/^^, = 1 - E'^/E^, and Y 

is the kinematic upper limit of y^c- The cross section 
asymmetry A^c for antiparallel (ti) and parallel (tt) ori- 
entations of the incoming muon and target electron spins 
is 



^'^^ - d^Ti 



dcrTi „ daTT 



dfjTT 



l-y^,/Y + y^,/2 
y^^l~y,./Y + ylJ2- ^''^^ 



The measured asymmetry Aoxp is related to A^o by 



^cxp(,y/jo j 



-fc -I /J ^pc yy^ic ) , 



(3.5) 

where P^ and P^ are the electron and muon polariza- 
tions, respectively. The measured asymmetries range 
from about 0.01 at low j/^o to 0.05 at high y^e- 



The experimental set-up for the fj,~e scattering mea- 
surement is shown schematically in Fig. ^(b). The lead 
foil is removed from the SVD and only the hodoscopes 
of the SVD are used to tag the incident muon which is 
tracked in three MWPC installed upstream of the mag- 
netized target. Between the target and the spectrometer 
magnet, three additional chambers measure the tracks of 
the scattered muon and of the knock-on electron. Down- 
stream of the magnet, the muon and the electron are 
tracked in two wire-chamber telescopes sharing a large 
MWPC. The electron is identified in the LGC and the 
muon is detected in a scintillation-counter hodoscopc lo- 
cated behind a 2 m thick iron absorber. 

The polarized electron target is a 2.7 mm thick foil 
made of a ferromagnetic alloy consisting of 49% Fe, 
49% Co and 2% V. It is installed in the gap of a soft- 
iron flat-magnet circuit with two magnetizing solenoidal 
coils iQ . The magnet circuit creates a saturated homo- 
geneous field of 2.3 T along the plane of the target foil. 
In order to obtain a component of electron polarization 
parallel to the beam, the target foil was positioned at an 
angle of 25° to the beam axis. 

To determine the target polarization, the magnetic flux 
in the foil under reversal of the target-field orientation is 
measured with a pick-up coil wound around the target. 
The magnetization of the target was found to be constant 
along the foil to within 0.3%. The electron polarization 
is determined from the magneto-mechanical ratio g' of 
the foil material. A measurement of g' for the alloy used 
does not exist; a value of g' — 1.916 ± 0.002 has been 
reported for an alloy of 50% Fe and 50% Co (7|]. We 
assume that the addition of 2% V does not affect g' but 
we enlarge the uncertainty to ±0.02. The resulting po- 
larization along the beam axis is \Pe\ = 0.0756 ± 0.0008. 
The loss of fj,-e events because of the internal motion of 
K-shell electrons Q affects the asymmetry A^xp by less 
than -0.001 and was therefore neglected. 

To measure the cross section asymmetry, the target- 
field orientation was changed between SPS pulses by 
reversing the current in the coil. The vertical compo- 
nent of the magnetizing field provides a bending power 
of 0.05 Tm which gives rise to a false asymmetry. This 
effect was compensated for by alternating the target an- 
gle every hour between 25° and —25° and averaging the 
asymmetries obtained with the two orientations. 

The trigger requires a coincidence between the two 
SVD hodoscope planes, an energy deposition of 15GeV 
or more in the LGC, and a signal in the muon ho- 
doscope (MH). The scattering vertex is reconstructed 
from the track upstream and the two tracks downstream 
of the magnetized target. The three tracks were required 
to be in the same plane to within 20° and the recon- 
structed vertex to be within ±50 cm of the target posi- 
tion. The two outgoing tracks were required to have an 
opening angle larger than 2 mrad and to satisfy the two- 
body kinematics of elastic scattering to within 1 mrad. 
Since the electron radiates in the target, we use the scat- 
tered muon energy to calculate y^o- 



11 



(U 

=1 



CO 




=1 



CO 



0.2 0.4 0.6 0.^ 




0.2 0.4 0.6 0.! 



FIG. 6. The QED radiative corrections to the asymmetry Afj,c (a) without experimental cuts, (b) The asymmetry if the 
following experimental cuts are included in the calculation: (i) recoil electron energy greater than 35 GeV, (ii) energy difference 
between initial and final states less than 40 GeV, (iii) angular cuts on both outgoing muon and electron. The corrections 5^c 
are given in percent. 



Background originates from brenisstrahlung (/^+ -^ 
/i+7) followed by conversion, and pair production (/i+ — > 
/i^e^e"). It was determined experimentally by using a 
/i~ beam with a similar set-up and triggering on /i~e~'' 
coincidences. Most of the background was eliminated by 
requiring that the energy conservation between the ini- 
tial and final states be satisfied within 40 GeV. This re- 
quirement rejects very few good events. The background 
correction to the beam polarization is —0.012 ± 0.004. 

The experimental asymmetry was obtained from data 
samples taken with the two different target field orien- 
tations. The data samples were normalized to the in- 
cident muon fiuxes using a random trigger technique. 
A possible false asymmetry due to the target magnetic 
field was studied using both a Monte Carlo simulation 
of the apparatus and data taken with an unpolarized 
polystyrene target under the same experimental condi- 
tions. In both cases the resulting asymmetry was found 
to be consistent with zero. The radiative corrections 
Sfj.e — ( -AQ™ /A^e — 1) to the first order cross section 
of Eq. ( |3.3D are evaluated using the program ^ela [ [73[ . 
The corrections are calculated up to C'(aQED) with finite 
muon mass and found to be negligible once the experi- 
mental cuts are apphed (Fig. ||). 

The polarization P^ = Acxpiy^J.c)/Af^ciy^lc) Pc in bins 
of y^e is shown in Fig. [^. The main contributions to 
the systematic error are the uncertainty of the fiux nor- 
malization, the false asymmetry, the uncertainty of the 
target polarization, and the background subtraction. 



3. The beam polarization 

The beam polarization obtained from the fi-e scatter- 
ing experiment in 1993 is |74,p3|: 



P., 



-0.779 ± 0.026 (stat.) ± 0.017 (syst.) (3.6) 



for E^ = 187.9 



GeV. The polarization measured by 
the muon decay method in 1993, P^ = —0.803 ± 
0.029 (stat.)±0.020 (syst.), has been published earlier [|. 
Both results are compatible. An alternative analysis 
with a larger data sample for the muon decay method 
is in progress and the systematic uncertainties of our 
previous analysis are be ing re-evaluated. The result of 
the /i-e scattering Eq. ( pig ) is used in this paper. For 



E.,. 



100 GeV a value of P^ = 



-0.82 ± 0.06 was used for 



the analysis of the A2 measurement. This is based on the 
measurement reported in Ref. |64|. Monte Carlo simu- 
lations of the muon beam [ pOf are consistent with these 
measurements of P^ for both beam energies. We have 
evaluated the average polarization of our accepted event 
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FIG. 7. Beam polarization vs. the ratio of electron and 
muon energies from polarized ^-e scattering. The dashed 
line represents the average value. 
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FIG. 8. Cross section of the SMC polarized target. 



sample taking into account the energy dependence of the 
muon polarization. The polarizati on w as calculated on 
an event-by-event basis using Eq. ( |3.l[ ) and assuming a 
monoenergetic pion beam (Fig. a). 



D. The polarized target 

The polarized proton target uses the method of dy- 
namic nuclear polarization (DNP) U& and contains two 
oppositely polarized target cells exposed to the same 
muon beam (Fig. H) B] . The solid target material is bu- 
tanol (CH3(CH2)30H) plus 5% water doped with para- 
magnetic EHBA-Cr(V) molecules. A superconducting 
magnet system |7^ and a •^He-^'^He dilution refrigera- 
tor (DR) [ra provide the strong magnetic field and the 
low temperature required for high polarization, and allow 
for frequent inversion of the field and thus of the polar- 
ization vectors. Additional subsystems include a double 
microwave set-up needed for the DNP and a 10-channel 
NMR system to measure the spin polarization |^ . Dur- 
ing data-taking, the nuclear spin axis is aligned either 
along or perpendicular to the beam direction in order to 
measure A\\ or A±, respectively. 

The two target cells were each 60 cm long, cylindrical. 



polyester-epoxy mesh cartridges of 5 cm diameter, sep- 
arated by a 30 cm gap. The target consisted of 1.8 mm 
butanol glass beads. The total amount of target ma- 
terial was 1.42 kg, with a packing fraction of 0.62 and 
a density of 0.985 g/cm^ at 77 K. The concentration of 
paramagnetic electron spins in the target material was 
6.2 X 10^^ spins/ml. In addition to butanol, the target 
cells contained other material, mostly the "^He-^He cool- 
ing liquid and the NMR coils for the polarization mea- 
surement (Table O . 

In the 2.5 T field and at a temperature below 1 K, the 
electron spins are nearly 100% polarized. When their 
resonance line is saturated at a frequency just above or 
below the absorption spectrum centered around the fre- 
quency of Ve ~ 69.3 GHz at 2.5 T, negative and positive 
proton polarizations are obtained. This technique was 
applied to polarize the material in the two target cells 
in opposite directions. Modulation of the microwave fre- 
quencies with a 30 MHz amplitude and a 1 kHz rate in- 
creased the polarization build-up rate by 20% and re- 
sulted in a gain in maximum polarization of 6%. This 
method was originally developed to improve the polar- 
ization of a deuterated butanol target Q • 

The DR |^l| cools the target material to a temperature 
below 0.5 K while absorbing the microwave power applied 
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for DNP. Once a high polarization is reached, the mi- 
crowaves are turned off and the target material is cooled 
to 50 mK. At this temperature the proton spin-lattice re- 
laxation time exceeds 1000 hours at 0.5 T. Under these 
'frozen spin' conditions, the polarization is preserved dur- 
ing field rotation and during measurements with trans- 
verse spin. To avoid possible systematic errors, the pro- 
ton polarizations were reversed by DNP once a week. 

The superconducting magnet system [|7^,Q consists of 
a solenoid with a longitudinal field of 2.5 T aligned with 
the beam axis, and a dipole providing a perpendicular 
'holding' field of 0.5 T. The solenoid has a bore of 26.5 cm 
into which the DR with the target cells is inserted; this 
diameter corresponds to an opening angle of ±65mrad 
with respect to the upstream end of the target. Sixteen 
correction coils allow the field to be adjusted to a relative 
homogeneity of ±3.5 x 10^^ over the target volume. In 
addition, the trim coils were used to suppress the super- 
radiance effect [Q , which can cause losses of the nega- 
tive proton polarization while the field is being changed. 
The spin directions were reversed every five hours with 
relative polarization losses of less than 0.2%. This was 
accomplished by rotating the magnetic field vector of the 
superimposed solenoid and dipole fields, with a loss of 
data-taking time of only 10 minutes per rotation Q. 
The dipole field was also used to hold the spin direction 
transverse to the beam for the measurement of A± . 

The proton polarization was measured with ten series- 
tuned Q-meter circuits with five NMR coils in each target 
cell gJ,Q. The polarization is proportional to the in- 
tegrated NMR absorption signal which was determined 
from consecutively measured response functions of the 
circuit with and without the NMR signal. The latter 
was obtained by increasing the magnetic field, and thus 
shifting the proton NMR spectrum outside the integra- 
tion window. The calibration constant was obtained from 
a measurement of the thermal equilibrium (TE) signals 
at 1 K, where the polarization is known from the Curie 
law Pte = ta.nh{hvp/2kT) ~ 0.002553; T is the lat- 
tice temperature, k the Boltzmann constant, and Vp is 
the proton Larmor frequency. The accuracy of the TE 
calibration signal contributed to the polarization error 
by AP/P = 1.1% |f|. The NMR signals were measured 
every minute during data-taking. The polarizations mea- 
sured with the individual coils were averaged for each tar- 
get cell and over the duration of one data taking run of 
typically 30 minutes. All measurements inside the same 
cell agreed to better than 3%. To detect a possible radial 
inhomogeneity, two of the five coils in the upstream tar- 
get cell were at the same longitudinal position, but one 
was in the center and the other at a radius of 1 cm. No 
significant difference was found between the polarizations 
measured by these two coils. 

The characteristic polarization build-up time was two 
to three hours. However, the highest polarizations of 
-1-0.93 and —0.94 were achieved only after several days 
of DNP. The average polarization during the data-taking 
was 0.86, and the relative error in the average polariza- 



TABLE II. Quantities (in moles) of the various chemical 
elements in the target volume. 



Element 


Quantity 


Element 


Quantity 


Element 


Quantity 


^H 


185.70 


F 


0.24 


Cu 


00.36 


^He 


6.00 


Na 


0.17 





22.70 


^He 


23.00 


Cr 


0.17 


c 


71.80 


Ni 


0.14 











tion of the target was estimated to be 3%. 



E. Muon spectrometer and event reconstruction 

The spectrometer is similar to the set-ups used by the 
EMC |86| and the NMC (Fig. ||). Aging chambers were 
replaced and new ones added to improve the redundancy 
of the muon tracking and to extend the kinematic cov- 
erage to smaller x. A major new streamer tube detector 
ST67 was constructed to identify and measure scattered 
muon positions downstream of the absorber. Triggers 
were optimized for improved kinematic coverage, in par- 
ticular in the region of small x. 



1. Spectrometer layout 

Three stages of the spectrometer can be distinguished: 
tracking of the incident muon, tracking and momentum 
measurement of the scattered muon, and muon identifi- 
cation. The beam tracking section upstream of the target 
is composed of two scintillator hodoscopcs (BHA/BHB) 
and the POB MWPC. A set of veto counters (VI. 5, V3, 
V2.1 and V2) defines the beam spot size. Beam tracks 
are reconstructed with an angular resolution of 0.1 mrad 
and an efficiency better than 90% for intensities up to 
5 X 10^ Ai/spill- 

The momentum of the scattered muon is measured 
with a conventional large-aperture dipole magnet (FSM) 
and a system of more than 100 planes of MWPC (Ta- 
ble III). The FSM is operated with bending powers of 
2.3 and 4.4 Tm at 100 GeV and 190 GeV beam energies, 
respectively, corresponding to a horizontal beam defiec- 
tion of 7 mrad. The angular resolution for scattered 
muons is 0.4 mrad. The large MWPC are complemented 
by smaller MWPC with a smaller wire pitch, to increase 
the redundancy and the resolution of the spectrometer in 
the high-rate environment at small scattering angles. 

Scattered muons are identified by the observation of 
a track behind a 2 m thick iron absorber. The muon 
identification system consists of streamer tubes, MWPC 
and drift tubes. To cope with the high beam intensity, 
the streamer tubes were operated with voltages at which 
their pulse heights were close to the electronic threshold. 
Their efficiencies were thus very sensitive to the ambient 
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TABLE III. Detectors of the muon spectrometer. 



Hodo- 


Modules 


Pitch 


Size 


Wire- 


Modules 


Pitch 


Size 


Dead 


scope 


X Planes 


(cm) 


(cm) 


chamber 


X Planes 


(cm) 


(cm) 


zone(cm) 


BHA-B 


2x8 


0.4 


8x8 


POA-E 


5x8 


0.1 


14 


— 


V123 


5x1 


— 


various 


PVl 


1x4 


0.2 


150x94 


— 


HI 


2 


7.0 


250x130 


PV2 


1x6 


0.2 


154x100 


8 


H2cal 


4 


28.0 


560x280 


P123 


3x3 


0.2 


180x80 


13 


H3 


2 


15.0 


750x340 


W12 


2x8 


2.0 


220x120 


12 


H4 


1 


15.0 


996x435 


W45 


6x4 


4.0 


530x260 


13-25 


Hl',3',4' 


1 


1.4 


50x50 


P45 


5x2 


0.2 


90 


12 


Sl,2,4 


1 


— 


various 


ST67 


4x8 


1.0 


410x410 


16 


H5 


1x2 


various 


19x20 


P67 


4x2 


0.2 


90 


12 


H6 


1x2 


various 


14 


DT67 


3x4 


5.2 


500x420 


83x83 



pressure and temperature, and a high-voltage feedback 
system was developed to stabilize the average streamer 
pulse height within 1%. 



2. Triggers 

The read-out of the detectors was triggered by prede- 
fined coincidence patterns of hits in different planes of 
scintillation-counter hodoscopes. Three physics triggers 
provide a coverage of different x and Q^ ranges (Fig. M). 
All triggers require that there is no hit in any of the 
beam-defining veto counters. 

The large-angle trigger Tl requires a coincidence pat- 
tern of the hodoscopes HI, H3 and H4. This trigger 
has a good acceptance for scattering angles 6 larger than 
20 mrad. Target pointing of the scattered muon is also 
required. The acceptance decreases for smaller angles, 
but extends to 6* « 3 mrad. The small-angle trigger T2 
uses the smaller hodoscopes HI', H3' and H4'. This trig- 
ger covers the range 5 mrad < 9 < 15 mrad. It has a 
more limited x-range than Tl. However, at a given x, 
T2 selects events with lower Q^ than Tl. A small- a; trig- 
ger T14 is provided by the SI, S2 and S4 counters which 
are placed close to the beam to cover scattering angles 
down to 3 mrad with good efficiency. The counters for T2 
and T14 were located on the bending side of the spec- 
trometer magnet. The acceptance of the triggers Tl and 
T14 extends down to x ~ 0.5 x 10"'^ and thus is sensi- 
tive to elastic scattering of muons from atomic electrons, 
X = rric/mp (Fig. H). The trigger rate per SPS spill was 
about 200 for Tl, 50 for T2 and 100 for T14. 

Other triggers include normalization and beam-halo 
triggers which were used for calibration, alignment, and 
efficiency calculations. 



3. Event reconstruction 

The track finding starts with the beam-track recon- 
struction. The nromentum of the incident muons is conr- 
puted from the hit pattern in the BMS hodoscopes. The 
beam track upstream of the target is found from the 
hits in the BHA and BHB hodoscopes and the FOB wire 
chamber. A coincidence is required between the hits in 
the BMS and those in the beam hodoscopes. 

The reconstruction of the scattered muon tracks starts 
in the muon identification system behind the hadron ab- 
sorber (ST67, DT67, P67). Tracks formd in this sys- 
tem are extrapolated upstream and reconstructed in the 
MWPC and drift chambers between the absorber and the 
FSM (W45, P45, W12, FOE). The next step in the recon- 
struction is the track finding in the FSM chambers (P123, 
FOD), starting with the vertical coordinates which are 
fitted by straight lines. Horizontal coordinates match- 
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FIG. 9. Kinematic ranges for triggers Tl, T2 and T14 at 
190 GeV. 
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TABLE IV. Kinematic cuts applied for the A\ 


and A± analysis. 




Kinematic 
variable 


j4|| analysis 
£^ = 190 GeV 






A± analysis 
£^ = 100 GeV 


V 

y 

e 


> 15 GeV 
< 0.9 

> 19 GeV 

> 9 mrad 








> 10 GeV 
< 0.9 

> 15 GeV 

> 13 mrad 




Final Data Sample for A\ 


analysis 


Final Data 


Sample for A± analysis 


X range 
Q^ range 

Events 


0.003 < X < 0.7 0.0008 < a: < 0.7 

1 < Q" < 90 0.2 < Q2 < 90 

4.5 X 10^ 6.0 X 10^ 


0.006 < 
1<Q2 

8.8 X 


x<0.6 
<30 
10^ 


0.0035 < X < 0.6 

0.5 < Q^ < 30 

9.6 X 10" 



ing the downstream tracks are searched for on circular 
trajectories inside the FSM. Because of the high track 
multiplicity in the FSM aperture, each extrapolation of 
a downstream track through the magnetic field is tested 
with a spline fit and the best track is retained. In the 
vertex chambers (PV12, POC), hits are selected using the 
extrapolated track reconstructed in the magnet, and are 
fitted by a straight line. It is verified that the recon- 
structed muon track satisfies the trigger conditions. 

The vertex position in the target is computed as the 
point of closest distance of approach between the beam 
and the scattered- muon tracks. Tracks are propagated 
through the magnetic field in the target using a Runge- 
Kutta method, taking into account energy loss and multi- 
ple scattering. In case of multiple beam tracks, the vertex 
with the best space-time correlation between the beam 
and the scattered-muon track is chosen. The vertex is 
reconstructed with resolutions of better than 30 mm and 
0.3 mm along and perpendicular to the beam direction, 
respectively. 



G. Event selection 

Since the A|| and A± data were recorded at different 
beam energies, they cover different kinematic ranges and 
are subject to different kinematic cuts (Table IV). A cut 
at small v rejects events with poor kinematic resolution, 
whereas a cut at high y removes events with large radia- 
tive corrections. A cut on the momentum of the outgoing 
muon reduces the contamination by muons from n and 
K production in the target and subsequent decay to a 
few 10"'^. The cut on was only applied for the analy- 
sis with Q^ >l GeV^. It rejects events with poor vertex 
resolution. 

Cuts were also applied to the beam phase space to en- 
sure that the beam fiux was the same for both target cells. 
Fiducial cuts on the target volume reject events from ma- 
terial outside the target cells (Fig. |l^). Less than 10% 
of the raw data were discarded because of instabilities in 
the beam intensity, detector efficiencies, and target po- 
larization. The size of the final data samples after all 
cuts is shown in Table IV. 



F. Data-taking 



IV. DATA ANALYSIS 



The data presented in this paper were taken during 
134 days of the 1993 CERN SPS fixed-target run. Most 
data were taken with longitudinal target polarization, at 
a beam energy of 190 GeV. For 22 days, data were taken 
with the target polarized transversely to the beam, at a 
beam energy of 100 GeV. 

A total of 1.6x10^ deep-inelastic-scattering events 
were reconstructed from the data with a longitudinally 
polarized target, using the three physics triggers Tl, T2 
and T14. The integrated muon flux was 1.7 x 10^"^. 

With transverse target polarization, only Tl was used 
and 1.6 million events were reconstructed. The trans- 
verse target field was always in the same vertical direc- 
tion and the spin direction was inverted by microwave 
reversal a total of 10 times. The integrated muon flux at 
100 GeV was 0.2 x lO^^. 



A. Evaluation of cross section asymmetries 



The two cross section asymmetries Au and Aj_ 



(Eq. (2/7)) are evaluated from counting rate asymmetries. 
To determine A\\ the four measured counting rates from 
the upstream and downstream target cells with the two 
possible antiparallel target spin configurations are used. 
The quantity At = A± cos 4> is determined separately for 
the upstream and downstream target cells from the four 
counting rates into the upper and lower vertical halves of 
the spectrometer for the two transverse spin directions. 



1. The yl|| analysis 

The number of muons N^ and A^d scattered in the up- 
stream and downstream target cells, respectively, is given 
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by 









(4.1) 
(4.2) 



fP^PdA^i), 

where $ is the integrated beam flux, P^ and Pd are the 
polarizations in the two target ceUs, riu and n^ the area 
densities of the target nucleons, and Cu and ad are the 
corresponding spectrometer acceptances. The dilution 
factor / accounts for the fact that on ly a f raction of the 
target nucleons is polarized (Section IV C ). The flux $ 
and the spin-independent cross section W cancel in the 
evaluation of the raw counting-rate asymmetries, Araw 
and ^J^AW obtained before and after target polarization 
reversal: 



iRAW 



^u + A^d ■ 



A' 



RAW 



K 



N' 



Provided that the ratio of acceptances is the same be- 
fore and after polarization reversal, i.e. a-a/a^ = a'^/a'^, 
and since n^/rid is constant, the acceptances a and the 
densities n cancel in the average of the raw asymmetries, 
so that 



A 



1 



fP^^Pt 



A 



RAW 



A' 



RAW 



If flu/ftd 7^ o.'u/o!d a 'false' asymmetry ensues, 



^falso — : 



1 



2fDP,,Pt 



r- 1 



r' -\ 
r' + l 



The 

photon-proton asymmetry A\ 

thus given by: 



(4.4) 



(4.5) 



virtual 



vlll/Li (Eq. ( |2.15D ) is 



Ai 



1 



IDPuPt 



^RAW + ^F 



RAW 



A 



false • 



(4.6) 



In these exp ressions, D is the depolarization fac- 
tor (Eq. ( 2.13 )), r — n^au/ndad, r' — nua[^/nda'^ and Pt 
is the weighted average of the target cell polarizations. 



2a 



ZN^ + ZNd 



EK + EK 



(4.7) 



Equation (4.6) provides an unbiased estimate of the cross 
section asymmetry for large numbers of events. To 
avoid possible biases for the number of events involved, a 
maximum likelihood technique was developed which al- 
lows a common analysis of all events in each x-bin. In 
this method, Au/D is computed from the event weights 
w = fDP^ using the expression 



, 1 






Ewd + E< 



A 



false ■ 



(4.8) 



As explained in Section IV C , in the actual analysis we 
use a weight w = f'DP^. A Monte Carlo simulation con- 
firmed that this method does not introduce any biases. 
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FIG. 10. Vertex distributions of scattered muons after 
kinematic cuts: (a) along the beam direction and (b) in the 
plane perpendicular to the target axis, at the location of one 
of the NMR coils. In (a), the dashed lines indicate the fidu- 
cial cuts on the target volume which coincide with the entry 
and exit windows of the target cells; most events outside the 
shaded region originate from interactions with the ^He-'*He 
cooling liquid. The small peak at a: ~ —3.9 m arises from 
scattering in the exit window of the target cryostat. In (b), 
the outer circle indicates the wall of the target cells, and the 
inner circle shows the radial cut applied. Scattering from the 
tubular NMR coils is clearly visible. 



2. The A±^ analysis 

A similar formalism applies to the measurement of 
the transverse asymmetry Aj^, where the event yields 
are given by N{(j)) = n'^aa{\ — /P^Pt cos(/)Aj_). Here, 
A± is obtained for each target cell separately from 

[V(0) - V(0- 7r)]/[V(0) + V(0- tt)] and A^/d becomes 



A^ 
d 



-I 



'^P^Pt] 

^false ; 



f J2fd cos ' 
\J2ifd cos 4 



f J2fdcos( 
\J2ifdcoa4 



(4.9) 
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where (Pt) is the average target polarization before and 
after reversal in absolute value. To obtain the same sta- 
tistical accuracy for A±/d and for A\\/D more data are 
required for A^/d due to its dependence on cos0, and 
also to a lesser extent to the fact that d < D. 



B. Radiative corrections 

QED radiative corrections are applied to convert the 

and (4.9) to one-photon ex- 



0.16 



measured asymmetries (4. 
change asymmetries, 
using: 



These corrections are calculated 



cr^ = va^'^ -I- CTtaii, 



(4.10) 



where a is the total, i.e. measured, spin-independent 
cross-section, a^'^ is the corresponding one-photon ex- 
change cross section, and o^taii is the contribution to ct 
from the elastic tail and the inelastic continuum. The 
corresponding differences of the cross sections for antipar- 
allel and parallel orientations of lepton and target spins 
are denoted by Act. The factor v accounts for vacuum 
polarization and also includes contributions from th e in- 
elastic tail close in x. The decomposition in Eq. 4.1C 
depends on the fraction of the inelastic tail included in v 
and is therefore to some extend ambiguous. Due to a can- 
celation of the different contributions, v is close to unity. 
Using the program TERAD ^ we find 0.98 <v < 1.03 
in the kinematic range of our data. For simplicity we set 
V to unity in our analysis and attribute all corrections to 

CTtaii [|3- 

Neglecting A2 and thus implying Ai = Ao'/(2I?cr), the 
radiative corrections to the one-photon asymmetry, A^ ^ 
can be written as 



Al 



p(^r+^r), 



(4.11) 



with p = wct^'^/ct'^ and Af = Aate,ii/2vDa'^'^. 

The ratio 'a^'^ j'a^ and the correction A" are evaluated 
using the program POLRAD ||89| , P0| . The asymmetry 
■^x{x) required as input is taken from Refs. [||JyJ§l and 
the contribution from J^ is neglected. The uncertainty in 
J^-^ is estimated by varying the input values of J^ within 
the errors. The factor p and the additive correction A'^ 
are shown in Table ^ at the average Q^ of each x-bin. 

We have incorporated p into the evaluation of the di- 
lution factor, /' = p/, on an event-by-event basis. Using 
the weight w — f'DPf^ we directly obtain Aj /p on the 
left-hand side of Eq. ^and thus AI'' (Eq. ([4.11D). 



The radiative corrections to the transverse asymmetry 
A"^ are evaluated as above, however assuming that 52 — 
92^ |p5| . The additive correction is much smaller than 
the statistical error and has been neglected. 




X 

FIG. 11. The dilution factor / (solid line) and the effective 
dilution factor /' = pf (dashed line) as a function of x. 



C. Dilution factor 

In addition to butanol, the target cells contain the 
NMR coils and the ^He -^He coolant mixture. The com- 
position in terms of chemical elements is summarized in 
Table y. The dilution factor / can be expressed in terms 
of the number ua of nuclei with mass number A and the 
corresponding total spin-independent cross sections a^ 
per nucleon for all the elements involved: 



,/ = 



J2a ^A ■ ^A 



(4.12) 



The total cross section ratios (j\/a^ for D, He, C 
and Ca are obtained from the structure function ratios 
^2V-Pf '^^^ Fi'/Fi @- The original procedure 
leading from the measured cross section rations '5\/'d% 
to the published structure function ratios was inverted 
step by step involving the isoscalarity corrections and 
radiative corrections (TERAD). For unmeasured nuclei 
the cross section ratios are obtained in the same way 
from a parameterization of ^2 (^)/^2'(^) ^^ ^ function 
of A ||-||. 

The dilution factor also accounts for the contamina- 
tion from material outside the finite target cells due to 
vertex resolution. This correction is applied as a function 
of the scattering angle, and the largest contamination oc- 
curs for the angles between 2 and 9 mrad, which results 
in a reduction of the dilution factor by about 6%. The 
correction needed because of the NMR coils (Fig. |l^) is 
convoluted with the distribution of the beam intensity 
profile. 



In the actual evaluation of Eqs. ( 4.8|) and (4.9) we use 
an effective dilution factor /' (Fig. 



/' = pf. 



(4.13) 
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TABLE V. The virtual photon-proton asymmetry A^ for Q^ > 1 GeV^ (top) and Q^ > 0.2 GeV^ (bottom). In the last 
colu mn, the first error is statistical and the second is systematic. (Araw) is the straight aver age of Araw and Araw in 



Eq. (4.4). The values for A^ have been corrected for radiative effects as described in Section 



[VB 



X range 



(x) 



(GeV") 



(P.) 



(y) 



(D) 



if) 



(p) 



(A-fiAw) 



A\'^ 



Al 



.003-.006 


.005 


1.320 


-.79 


.791 


.80 


.070 


1.50 


.004 


.007 


.083±.041±.006 


.006-.010 


.008 


2.068 


-.78 


.748 


.76 


.081 


1.39 


.003 


.008 


.044±.037±.004 


.010-.020 


.014 


3.562 


-.78 


.704 


.72 


.090 


1.30 


.003 


.010 


.061±.032±.004 


.020-.030 


.025 


5.733 


-.78 


.660 


.68 


.096 


1.24 


.003 


.012 


.068±.044±.005 


.030-.040 


.035 


7.797 


-.78 


.634 


.66 


.099 


1.21 


.002 


.015 


.041±.052±.003 


.040-.060 


.049 


10.445 


-.78 


.603 


.64 


.102 


1.18 


.006 


.017 


.104±.045±.007 


.060-.100 


.077 


15.011 


-.78 


.551 


.60 


.106 


1.14 


.009 


.020 


.180±.045±.013 


.100-.150 


.122 


21.411 


-.78 


.498 


.55 


.112 


1.10 


.013 


.022 


.289±.058±.019 


.150-.200 


.173 


27.799 


-.79 


.456 


.51 


.118 


1.08 


.012 


.022 


.276±.080±.019 


.200-.300 


.242 


35.542 


-.79 


.417 


.47 


.127 


1.05 


.010 


.019 


.246±.082±.017 


.300-.400 


.342 


45.453 


-.78 


.377 


.43 


.139 


1.02 


.021 


.010 


.499±.132±.036 


.400-.700 


.482 


57.089 


-.78 


.337 


.39 


.156 


0.99 


.022 


-.006 


.527±.174±.041 


.0008-.0012 


.001 


0.285 


-.78 


.808 


.85 


.044 


1.74 


-.001 


.002 


-.032±.077±.004 


.0012-.002 


.002 


0.445 


-.78 


.794 


.83 


.054 


1.65 


.002 


.003 


.085 ±.055±.007 


.002-.003 


.003 


0.686 


-.78 


.781 


.80 


.062 


1.56 


.001 


.004 


.031 ±.054±.004 


.003-.006 


.004 


1.193 


-.78 


.763 


.77 


.073 


1.46 


.003 


.006 


.059 ±.034±.005 


.006-.010 


.008 


2.038 


-.78 


.738 


.75 


.082 


1.38 


.003 


.008 


.050 ±.036±.004 



as discussed in Section [VB| . The present procedure guar- 
antees a proper calculation of the statistical error in the 
asymmetry, in contrast to our previous analysis |9|-|lJ] 
where all radiative effects were included as an additive 
radiative correction. We find an increase in the statis- 
tical error by a factor 1/p which reaches 5 at small- a; 
(Table 0). However, the central values of the asymme- 
tries remain unaffected by the change in the radiative 
correction procedure [ p7[ . 

The dilution factor is shown in Fig. O where it is com- 
pared to the 'naive' expectation for a mixture of 62% 
butanol, (CH3(CH2)3 0H), and 38% helium by volume, 
/ ~ 0.123. The rise of / at x > 0.3 is due to the de- 
crease of the ratio F2/F2, whereas the drop in the low 
x-range is due to the larger contribution of radiative pro- 
cesses from elements with mass number much larger than 
hydrogen. 



D. The longitudinal cross section asymmetry 

1. Results for A^ 

The virtual photon asymmetry A^ is calculated from 
|4.8|), (4.11) and (4.13) under the assumption that 



Eqs 

^faiso — 0- The uncertainty int roduced by this assump- 
tion is estimated using Eq. (4.5). 

The results for J^fov Q^ > 1 GeV^ are shown in Ta- 




and in Fig. |l^. The kinematic quantities in Ta- 
are mean values within the bins calculated with 



the weighting factor (f'DP^). In addition to the re- 



sults given in Rcf. |9|, we in clude h ere data obtained 
with the T14 trigger (Section |IIIE2| ). In Table ^ and 
in Fig. |l3, we also show data in the kinematic range 
0.2 GeV^< Q2 < 1 GeV2, 0.0008 <x< 0.003. These 
data are not used to evaluate g^ or its first moment. 

The sources of systematic errors in A^ are time- 
dependence instabilities of the acceptance ratios r and 
r', uncertainties in the beam and target polarizations, in 
the effective dilution factor /', the radiative corrections, 
and in i? = g^/gy, and the neglect of A2 . The individual 
errors (Table VI) are combined in quadrature to obtain 






0.8 
0.6 
0.4 
0.2 

-0.2 



• Q^> I GeV^ 






Q^> 0.2 GeV^ 
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T 



10 



10 
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X 



FIG. 12. The virtual photon asymmetry A'^ as a function 
of X. The error bars show statistical errors only; the system- 
atic errors are indicated by the shaded area. 
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FIG. 13. The virtual photon-proton asymmetry A^ as a 
function of Q^, for constant values of x. The closed circles 
are data from this experiment. The data of the EMC and 
E143 experiments are also shown as open circles and squares, 
respectively. 



TABLE VI. Contributions to the systematic errors at the 
average Q^ of the 2;-bin. 



(x) A^M.e APt AP^ A/' Arc AA2 



AR 



0.005 


0.0021 


0.0025 


0.0033 


0.0016 


0.0012 


0.0006 


0.0027 


0.008 


0.0019 


0.0013 


0.0017 


0.0008 


0.0012 


0.0007 


0.0012 


0.014 


0.0019 


0.0018 


0.0024 


0.0011 


0.0011 


0.0009 


0.0021 


0.025 


0.0018 


0.0020 


0.0027 


0.0013 


0.0010 


0.0002 


0.0031 


0.035 


0.0018 


0.0012 


0.0016 


0.0008 


0.0010 


0.0003 


0.0016 


0.049 


0.0018 


0.0031 


0.0041 


0.0020 


0.0009 


0.0003 


0.0040 


0.077 


0.0019 


0.0054 


0.0071 


0.0035 


0.0009 


0.0004 


0.0080 


0.122 


0.0019 


0.0087 


0.0114 


0.0058 


0.0010 


0.0005 


0.0112 


0.173 


0.0020 


0.0083 


0.0109 


0.0056 


0.0010 


0.0005 


0.0110 


0.242 


0.0020 


0.0074 


0.0097 


0.0051 


0.0009 


0.0022 


0.0105 


0.342 


0.0020 


0.0150 


0.0197 


0.0107 


0.0007 


0.0025 


0.0236 


0.482 


0.0020 


0.0158 


0.0208 


0.0117 


0.0008 


0.0030 


0.0293 


0.0011 


0.0032 


0.0010 


0.0013 


0.0011 


0.0009 


0.0005 


0.0017 


0.0016 


0.0027 


0.0025 


0.0034 


0.0026 


0.0010 


0.0008 


0.0035 


0.0025 


0.0024 


0.0009 


0.0012 


0.0009 


0.0011 


0.0011 


0.0012 


0.0044 


0.0021 


0.0018 


0.0024 


0.0014 


0.0012 


0.0007 


0.0023 


0.0078 


0.0020 


0.0015 


0.0020 


0.0010 


0.0012 


0.0008 


0.0015 



2. Comparison with earlier experiments 

In Figure l3, we compare our results for A^ with data 
from earlier experiments |l|J^,|6|,^. Good agreement is 
observed in the kinematic region of overlap. A consis- 
tency test between the SLAG E80/E130, EMC, SLAG 
E143 and SMG data yields a x^ = 11.4 for 16 degrees of 
freedom. Since the average Q^ of SMG and E143 differ by 
a factor of seven, the good agreement confirms the earlier 
conclusion that no Q^ dependence is observed within the 
present accuracy of the data. 



E. The transverse cross section asymmetry 



the total sys tematic error (Table ^). 

Table ^ and Fig. |l| show A^ as a function of Q^ 
and X, including the data with Q^ < 1 GeV^. In Fig- 
a small correction is applied to the data to display 
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ure 

them at the same average x in each bin. A study of the 
Q^ dependence which includes the SMG data HJl^ was 
first made by the E143 collaboration for 0.03 < x < 0.6 
and Q^ > 0.3 GeV^, and showed no significant Q^ de- 
pendence for Q2 > 1 GeV^ |||]. We study here the Q^ 
dependence for 0.003 < x < 0.03. A parametrization 
Ai = a + blogQ^ is fitted to the data and b is found 
to be consistent with zero for all x in this range. When 
fitting a parametrization a' + c/Q^ to account for possi- 
ble higher twist effects, we again find no significant Q^ 
dependence. 



1. Results for A^ 

The asymmetry A^ is obtained from o ur me asurements 
of A\ and of ^p |i|,|,|, using Eq. ( |2l7| ). It is seen 

from Eq. (2.E), that A2 has an explicit l/vQ^ depen- 
dence and hence it is convenient to evaluate ^JQ^A^ as- 
suming that it is independent of Q^ in Eq. ( |4.9| ). Our 
results do not depend on this assumption [ p7[ . 

The re sults for the asymmetry A^ are shown in Ta- 
ble VIII and in Fig. na. They are significantly smaller 
than the positivity limit \A2\ < yR and are consistent 
with A^ = and with the assumption that 172 = 92 ^ i 
i.e. 92 = 0. Also shown in Fig. |l| are the E143 data ^^. 
They confirm our results, with better statistical accuracy, 
for x > 0.03. 

The main systematic uncertainties are due to the 
parametrizations of A?./D and R. The effects due to time 
variations of the acceptance are negligible as expected. 
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TABLE VII. The virtual photon-proton asymmetry A^ as 
a function of x and Q^. Only statistical errors are shown. 



(x) 



(GeV^) 



A1 



(x) 



(GeV^) 



Al 



0.0009 
0.0010 
0.0011 



0.25 
0.30 
0.34 



-0.122± 
0.033 ± 
0.082 ± 



0.110 
0.137 
0.169 



0.0345 
0.0359 



7.77 
10.15 



0.058 ± 
~0.012± 



0.082 
0.095 



0.0014 


0.38 


0.209 ± 0.081 


0.0017 


0.46 


0.042 ± 0.102 


0.0018 


0.55 


-0.086± 0.109 


0.0023 


0.58 


0.114 ± 0.085 


0.0025 


0.70 


-0.009± 0.094 


0.0028 


0.82 


-0.025± 0.102 



0.0474 
0.0473 
0.0478 
0.0484 
0.0527 



2.94 
5.49 
7.83 
10.96 
14.73 



-1.114± 
~0.117± 
0.241 ± 
0.123 ± 
0.058 ± 



0.589 
0.142 
0.094 
0.068 
0.098 



0.0036 
0.0043 
0.0051 
0.0057 



0.88 
1.14 
1.43 
1.70 



-0.006± 
0.089 ± 
0.119 ± 

-0.033± 



0.065 
0.054 
0.067 
0.118 



0.0738 
0.0744 
0.0751 
0.0762 
0.0855 



5.33 
7.88 
11.09 
16.32 
23.04 



0.359 ± 
0.212 ± 
0.214 ± 
0.203 ± 
0.066 ± 



0.239 
0.142 
0.088 
0.068 
0.105 



0.0070 
0.0072 
0.0077 
0.0085 
0.0092 



1.42 
1.76 
2.04 
2.34 
2.72 



0.037 ± 
0.014 ± 
-0.045± 
0.166 ± 
0.145 ± 



0.094 
0.073 
0.071 
0.085 
0.093 



0.1193 
0.1199 
0.1204 
0.1208 
0.1293 



7.36 
11.16 
16.47 
24.84 
34.28 



0.456 ± 
0.480 ± 
0.364 ± 
0.199 ± 
0.172 ± 



0.242 
0.159 
0.110 
0.098 
0.137 



0.0122 
0.0125 
0.0141 
0.0165 
0.0184 



2.15 
2.82 
3.52 
4.43 
5.43 



0.184 ± 
0.020 ± 
0.066 ± 
0.085 ± 
-0.042± 



0.090 
0.067 
0.053 
0.069 
0.113 



0.1713 


14.15 


0.288 ± 0.143 


0.1717 


24.92 


0.349 ± 0.156 


0.1742 


39.54 


0.212 ± 0.123 



0.2384 
0.2396 
0.2462 



14.53 
29.71 
52.76 



0.139 ± 
0.110 ± 
0.413 ± 



0.176 
0.132 
0.131 



0.0235 
0.0236 
0.0242 
0.0263 



2.95 
4.38 
5.75 
7.42 



0.189 ± 

-0.026± 

0.107 ± 

0.072 ± 



0.176 
0.086 
0.070 
0.080 



0.3392 
0.3408 
0.3432 



15.29 
29.82 
61.49 



0.644 ± 
0.814 ± 
0.333 ± 



0.354 
0.241 
0.179 



0.0339 
0.0341 



4.14 
5.81 



0.003 ± 
0.097 ± 



0.174 
0.119 



0.4747 
0.4858 



26.74 
71.58 



0.541 ± 
0.518 ± 



0.306 
0.213 



TABLE VIII. Results on the asymmetry Ag. Only statis- 
tical errors are given. The A2 values are the average values 
from the two target cells. 



X range 



(x) 



(Q^) (GeV^ 



Al 



0.006 - 


- 0.015 


0.010 


1.4 


0.002 ±0.109 


0.015 - 


- 0.050 


0.026 


2.7 


0.041 ± 0.076 


0.050 - 


- 0.150 


0.080 


5.8 


0.017 ±0.099 


0.150 - 


- 0.600 


0.226 


11.8 


0.149 ±0.161 



0.0035 - 0.006 
0.006 - 0.015 



0.005 
0.01 



0.7 
1.3 



-0.066 ±0.167 
0.086 ± 0.097 



a, ^ 

^0 



1.0 
0.8 
0.6 
4 
0.2 
0.0 



■ E80/E130 
A EMC 

• SMC: Q^ >lGeV2 

o E143 : 29 GeV 
fl, 16 GeV 
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FIG. 14. The virtual photon-proton asymmetry A'^ as a 
function of x from this experiment, compared with data from 
the EMC and the SLAG E80, E130, and E143 experiments. 
For E143, the structure function ratio g\/Ff is shown instead 
of A\. The errors are statistical only. 
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FIG. 15. Results for the asymmetry ^2(^) extrapolated to 
Qo — 5GeV^ assuming y'Q^^2 scales [hol. The solid and 
dashed curves show the limit \A2\ < VR and the prediction 
corresponding to 172 = 0, respectively. Also shown are data 
from the E143 experiment mn extrapolated to the same Qg 
assuming that \/cpA2 scales. The errors are statistical only. 



since the results depend on the ratio of acceptances for 
niuons scattered into the top and the bottom halves of 
the spectrometer, which should be affected in the same 
way by typical variations of chamber efficiencies. The 
errors from the dilution factor and the beam and target 
polarizations are also very small. The total systematic 
error on A^ is at least one order of magnitude smaller 
than the statistical error at all values of x. 



V. RESULTS FOR gl AND ITS FIRST MOMENT 

A. Evaluation oi g\{x,Q ) 

The spin-dependent structure function g^ is evalu- 
ated from the virtual photon-proton asymmetry A^ us- 
ing Eqs.( p.l5[) and (2.16). This analysis is restricted to 
Q^ > 1 GeV^ . For F2 , we use the parametrization of 
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TABLE IX. Results for the spin-dependent structure function gf. The first error is statistical and the second is systematic. 
The third error in the last column is the uncertainty associated with the QCD evolution. 



x-range 


(x) 


(Q'iGeY^)} 


9n^,Q') 


g1{x,Ql = WGeY^) 


0.003-0.006 


0.005 


1.3 


1.97±0.97±0.15 


2.37±0.97±0.15±0.66 


0.006-0.010 


0.008 


2.1 


0.73±0.61±0.06 


1.03±0.61±0.06±0.17 


0.010-0.020 


0.014 


3.6 


0.63±0.33±0.05 


0.79±0.33±0.05±0.04 


0.020-0.030 


0.025 


5.7 


0.45±0.29±0.03 


0.51±0.29±0.03±0.02 


0.030-0.040 


0.035 


7.8 


0.20±0.26±0.02 


0.22±0.26±0.02±0.01 


0.040-0.060 


0.049 


10.4 


0.38±0.17±0.02 


0.37±0.17±0.02±0.00 


0.060-0.100 


0.077 


15.0 


0.42±0.10±0.02 


0.40±0.10±0.02±0.01 


0.100-0.150 


0.122 


21.4 


0.41±0.08±0.03 


0.39±0.08±0.02±0.01 


0.150-0.200 


0.173 


27.8 


0.26±0.08±0.02 


0.25±0.08±0.02±0.01 


0.200-0.300 


0.242 


35.5 


0.15±0.05±0.01 


0.15±0.05±0.01±0.01 


0.300-0.400 


0.342 


45.5 


0.15±0.04±0.01 


0.17±0.04±0.01±0.00 


0.400-0.700 


0.482 


57.1 


0.06±0.02±0.00 


0.08±0.02±0.00±0.00 



Ref. |9^ and for R the parametrization of Ref . |£^ . The 
parametrization of R is based on data for x > 0.01 only 
and therefore must be extrapolated to cover smaller val- 
ues of X. However, the structure function gi at the av- 
erage Q^ of the measurement is nearly independent of R 
due to a partial cancelation between the R dependence 
of D, of F2, and of the explicit term (1 -I- R{x, Q^)). The 
results for gi are shown in Table [X and, together with 
our deuteron data |13[, in Fig. 16. 



B. Evolution of gf to a fixed Qq 

To evaluate the first moment T^ = J„ (/fda;, the mea- 
sured gi{x, Q^) must be evolved to a common Qg for all 
X. In previous analyses, gi(x, Qq) was obtained assuming 
Ai ~ gi/Fi to be independent of Q^. This assumption 
is consistent with the data. However, perturbative QCD 
predicts the Q^ dependences of gi and Fi to differ by a 
considerable amount at small-x. The evolution oi gi/Fi 
is poorly constrained by the data in this region, where 
the data cover a very narrow Q^ range. Recent exper- 
imental and theoretical progress allows us to perform a 
QCD analysis of polarized structure functions in next- 
to-leading order NLO, and therefore a realistic evolu- 
tion of gi can be o btained. Three groups have published 
such analyses pi, 100, 101 1. They all use the splitting 
and coefficient functions calculated to NLO in the MS 
scheme p3|-p5[, but the choices made for the reference 
scales Q^ at which the polarized parton distributions 
are parametrized and the forms of the parametrization 
are different. Also the selections of data sets used for the 



fits differ. In ReL [gl 
tions are transformec 



the splitting and coefficient func- 
from the MS scheme to different 
factorization schemes before the fits are performed. We 
shall refer to the results obtained in the Adler-Bardeen 
scheme. 
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FIG. 16. The structure functions g\ and gf at the mea- 
sured Q^ and the corresponding gr". The upper and lower 
shaded areas represent the systematic error for g\ and g'^, 
respectively. 



We used the method 
data and those of Refs. 



of RcL 1 31 1 to fit the present 
l|-|l|,|,|6f,§. The quark sin- 
glet, non-singlet and gluon polarized distributions are 
parametrized as 

A/(x, g?rf) = NfT^fx'^' (1 - xf' (1 + afx) , (5.1) 

where the normalization factors Nf are chosen such that 
/ A/dx = rjf. We have assumed that Og — oas- The 
normalizations of the non-singlet quark densities are fixed 
using neutron and hyperon (3 decay constants and assum- 
ing SU(3) flavor symmetry. We use gA/gv = F + D = 



^The code was kindly provided by the authors. 
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TABLE X. Parameters of the polarized parton distributions at Q^^i 
text. 



1 GeV , obtained from the QCD fit discussed in the 



Aq 

AE 



Tre~ 



25.4 ±39.1 


-0.67 ±0.25 


2.12 ±0.28 


proton: 1.087 ±0.006 (fixed) 
deuteron: 0.145 ±0.002 (fixed) 


1.30 ±0.16 


0.71 ±0.33 


1.56 ±1.00 


0.40 ± 0.04 


HAS 


-0.70 ±0.27 


4 (fixed) 


0.98 ±0.61 



-1.2601 ± 0.0025 [[02[ and F/D = 0.575 ± 0.016 p3 |. 

The parameters of the polarized parton distributions ob- , 

tained from this fit are given in Table ^ and the fit is 

shown in Fig. O. We have fixed the exponent /3 of the 

gluon dist ribution to /3 = 4 as expected from QCD count- 2 

ing rules |104|Jl05| ], while the fitted values of /3 for the q^^ 

quark singlet and non-singlet components are found to OjQ 

be close to the expectation /3 = 3. The x^ ^01 the fit is i 

284 for 295 degrees of freedom. It is important to note, 

however, that the fit does not converge without our data 

points for x < 0.03, where the Q^ range is narrow. Re- q 

suits of E142 on gf were not included in the fit, but used 

as a cross check. In Figure [l7| their data and g'^ ' calcu- o.4 

lated from the fit to g\ and g^ are presented, and found 
to be in very good agreement. 

The measured gi(x, Q^) are then evolved from Q^ to Oo 

Qq by adding the correction 



fit/ 



5g^{x, Q\ Qi) = gf{x, Qf,) - gf{x, Q') 



(5.2) 



where g^* is calculated by evolving the fitted parton dis- 
tributions. The resulting g\{x^ Q^) is shown in Table IX 
and Fig. 11^. Also show n is the 9^ (a:, Qq) obtained by us- 
ing the fits of Refs. |^, 100, 101 1, and by assuming scahng 
for gi/Fi. For the lowest x bin, the latter results in a 
considerably larger value of gi . 



C. The first moment of gf 



0), its first 



From the evolved structure function gi{x 
moment T^ is evaluated at Qq ~ 10 GeV^, which is close 
to the average Q^ of our data. The integral over the 
measured x-range is 

pO.7 

/ fff(a;, go)da; = 0.130 ± 0.013 ± 0.008 ± 0.005 , (5.3) 

io.003 

where the first error is statistical, the second systematic 
and the third is the uncertainty due to the Q^ evolution. 
The individual contributions to the systematic error are 
summarized in Table XI . The error from the evolution is 



mainly due to the uncertainties in the factorization and 
renormalization scales, in the parametrizations chosen for 
the parton distributions, the error in as{Mz) and mass 
threshold effects. In addition we varied the values of F 
and D used as inputs to the fit, and of the A2, Afaisc, 
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FIG. 17. The structure functions g\^ of, and g\ at the 
measured Q^ for the SMC ||l|, E143 |^ and E142 |07l 
data. The solid curves correspond to our NLO fits at the Q 
of the data points, the dashed curve at Q\=\{i GeV^, and the 
dot-dashed curve at Q\=\ GeV^. 



/, P^, Pt^ Fi and of the radiative corrections used to 
calculate g\. The uncertainty in the fitted parameters 
of the parton distributions is also included, but is found 
to be relatively small. These errors on bgiix^Q"^ are 
treated as correlated from bin to bin, but uncorrelated 
amongst each other. 

The resulting g\ using the different phenomenological 
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FIG. 18. The structure function g^ evolved to Qq — 10 
GeV^ using the scaling assumption that gi/Fi is independent 
of Q^, and using NLO evol ution accordin g to our analysis and 
those of BFR [ml, GRSV M and GS [101 . 



given in Eq. ( p.5[ ), while the x behavior in the diffractive 

model, gi{x) « (xln x)~^, gives /^ ^ gi{x,QQ)dx — 
0.036 ± 0.016. This model results in a larger T^, but 
cannot simultaneously accommodate the negative values 
of g" found from our combined deuteron ||l3] and pro- 
ton data (Fig. |l^). In principle the low- a; contribution 
to the integral can be obtained from the fit to gi, i.e. 
gf^. However, as known from unpolarized parton distri- 
bution functions, the behavior of the fitted distribution 
below the measured region is unreliable since it depends 
strongly on the choice of the function, renormalization, 
and factorization scales. 

The result for the first moment of gi{x, Qq) is 

rP(Q2 ^ l0GeV2) = 0.136±0.013±0.009±0.005 . (5.6) 

Using the resu lts of the NLO evolutions of Refs. [ pl| , 
| iOO| and |Oltl we find r{{Q^) between 0.133 to 0.136 



(Fig. 



If we evaluate 5^(2;, Qq) assuming that 51/-F1 



analyses of the Q^ evolution pi, IOC, 101 1 are shown in 
Fig. Hq. Despite their different procedures, the differences 
in their results are small and are covered by the error that 
we quote for the evolution uncertainty. 

To estimate the integral for 0.7 < x < 1.0 we assume 
that A^ = 0.7±0.3 in this region. This is consistent with 
the high-x data and with the expectation from perturba- 
tive QCD that 51/Fi ^ 1 as x ^ 1 pOl. We obtain 



1.0 



is independent of Q^ we obtain Tl{Ql) = 0.139±0.014± 
0.010. We conclude that within the experimental accu- 
racy of our data the different NLO QCD analyses yield 
consistent results for the evolution oi gi, and that gi/Fi 
deviates significantly from scaling at small x. 



D. Combined analysis of r!j' 

The combined analysis of Fj' includes the proton spin 
asymmetries for Q^ > 1 GeV^ from our data and those 
of' Refs. 0,1,1 shown in Fig. |l4l The EMC and SMC 
data were taken at an average Q^ of 10 GeV^, while for 
the SLAC data the average Q^ is 3 GeV^ . The combined 



0.7 



glix, Qt = 10 GeV^)dx = 0.0015 ± 0.0007 . (5.4) ^^^^^j, ^^ ^^^j^^,^^ ^, an intermediate Q' of 5 GeV^ to 



The results from our fit shown in Fig. |l7| are used to 
evaluate L „„., .g^(a;, Q^)da; and found to be consistent 
with the sum of Eqs. (^^ and ( [s!^ ). 

The contribution to the first monrent from the unmea- 
sured region < x < 0.003 is evaluated assuming a con- 
stant g^ a,t Q^ — 10 GeV^, in agreement with a Regge- 
type behavior [g^ • Using the average of the two lowest 
X data points in Table [X we obtain 



0.003 



gP(x, Ql = 10 GeV2)da; = 0.0042 ± 0.0016 . (5.5) 



However, to evaluate the systematic error on F^* we have 
assumed an error of 100% in this integral (Table XI). It 
should be noted that we have assumed constant Regge- 
type behavior at Q^ — lOGeV^. If we apply the same 
procedure at Q^ = 1 GeV^ and then evolve the result- 
ing extrapolation to Q^ = 10 GeV^ using the NLO fits, 
we obtain a value which is within 1.5(7 of the assumed 
error. Oth er m odels describing the small- a; behavior of 
gi (Section p^ID ) are also considered to check the sen- 
sitivity of our result to the small- a; extrapolation. A 
gi{x) « Ina; dependence is compatible with the error 



avoid a large Q^ evolutions. Corrections to gi/Fi cal- 
culated at NLO are found to be up to 20-25%. The 

TABLE XL Contributions to the error of V^ 

Source of the error AFi 

Beam polarization 

Extrapolation at low x 

Target polarization 

Uncertainty on F2 

Dilution factor 

Acceptance variation Ar 

Momentum measurement 

Kinematic resolution 

Radiative corrections 

Extrapolation at high x 

Neglect of A2 

Uncertainty on R 



0.0048 
0.0042 
0.0036 
0.0030 
0.0025 
0.0014 
0.0014 
0.0010 
0.0008 
0.0007 
0.0004 
0.0000 



Total systematic error 


0.0087 


Evolution error 


0.0045 


Statistical error 


0.0125 
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TABLE XIL r^" and the contributions from different x regions at Qg =5 GeV^. The results of our analysis of the SMC and 
the E143 data, as well as the combined analysis of the SLAC-E80/130 0, EMC ||], SMC and SLAC-E143 § data are given 
with the statistical and systematic errors added in quadrature. Results of extrapolations are marked with an (*). 



X range 


0-0.003 


0.003-0.03 


0.03-0.7 


0.7-0.8 


0.8-1 


0-1 


SMC 
EMS 
ALL 


0.004(2)* 
0.0012(1)* 
0.004(2)* 


0.022(7) 

0.010(1)* 

0.021(6) 


0.104(13) 
0.115(7) 
0.114(6) 


0.0018(4)* 

0.0020(6) 

0.0020(6) 


0.0006(2)* 
0.0006(2)* 
0.0006(2)* 


0.132(17) 
0.129(8) 
0.141(11) 



TABLE XIIL The EUis-Jaffe sum rule calculated with NLO QCD corrections compared to our result for F^ at Q§ =10 GeV^ 
and 5 GeV^ and to the combined analysis of fit the E80/E1S0 §, EMC (|, SMC and E143 § data at Ql =5 GeV^ The 
Bjorken sum rule calculated with NNLO QCD corrections and compared to our results on V\ — FJ from the SMC, the combined 
analysis of V\ and Ff (SMC ||l| and EMS §) and the combined analysis of Ff ,F5' and F? (E142 fff^). 



Exp eriment /Theory 


VI 


ri 


r? 


V\ - F? 


Qt = 10 GeV" 


SMC 
EUis-Jaffe/Bjorken 


0.136 ±0.016 
0.170 ±0.004 


-0.046 ± 0.021 
-0.016 ±0.004 


0.041 ±0.007 
0.071 ±0.004 


0.183 ±0.034 
0.187 ±0.002 


Qg = 5 GeV" 



SMC 

COMBINED(p,d) 

COMBINED(p,d,n) 

EUis-Jaffe/Bjorken 



0.132 ±0.017 
0.141 ±0.011 
0.142 ±0.011 
0.167 ±0.005 



-0.048 ± 0.022 
-0.065 ±0.017 
-0.061 ±0.016 
-0.015 ±0.004 



0.039 ± 0.008 
0.039 ± 0.006 
0.038 ± 0.006 
0.070 ± 0.004 



0.181 ±0.035 
0.199 ±0.025 
0.202 ± 0.022 
0.181 ±0.003 



evolution of g\ to Q^ — 5 Gc V^ (F ig. 19) is performed 



The combined result for the first moment of g\ is 



using the procedure of Section VB. 

The data are combined on a bin- by- bin basis. The in- 
tegrals AF^^ = J^^ gi{x; Q^) dx are computed for the x- 
bins of each experiment individually, starting from the 
published asymmetries. The AF^^ which fall into the 
same SMC x-bin are first summed for each experiment 
and then the integral for this bin is obtained as weighted 
average of these sums. The weights are calculated by 
adding the statistical errors and systematic errors un- 
correlated between the experiments in quadrature. The 
error and the central value of the integral in the mea- 
sured region is computed using a Monte Carlo method, 
which takes into account the bin-to-bin correlation of the 
systematic errors within each experiment as well as cor- 
relations between the experiments. These correlated con- 
tributions are due to the polarizations of the beam and 
the target, the dilution factor, the neglect of ^2, the time 
dependence of the acceptance ratio, the radiative correc- 
tions, and the parametrizations of F2 IQ , of R |9^ , and 
of the parton distribution functions used to evolve gi. 
Correlations between the experiments arise mainly from 
the latter three sources. The error distributions in the 
Monte Carlo sampling are assumed to be Gaussian. 

The X range of the combined data is 0.003 < a: < 0.8. 
The extrapolations at large and small x are performed 
using the procedures described in Section V C . The con- 
tributions to the integral from the mea sure d and extrap- 
olated regions of x are shown in Table XII . 



T\{Ql = 5 GeV^) = 0.141 ± 0.011 



{All proton data) . 
(5.7) 



If Ai is assumed to be independent of Q^, we obtain 
T\ = 0.140 ±0.012. 

It should be noted that the error quoted by the E143 
collaboration p] from their data alone and the error ob- 
tained from our combined analysis are comparable. The 
statistical uncertainties of the SMC data for 0.003 < a; < 
0.03 introduce a larger error to F^ than the uncertainty 
assumed by the E143 collaboration for their extrapola- 
tion from X = 0.03 to x = 0. We also calculated the 
extrapolations from the evolved E143 and SMC data sep- 



arately. The results are compared in Table XII 



The results for F^ from SMC and from the combined 
analy sis ar e compared with the EUis-Jaffe sum rule in 
Table XIII . The EUis-Jaffe prediction is calculated from 
Eq. (2.44). The higher-order QCD corrections are ap- 
plied assuming three active quark flavors, and using 
as{5 GeV^) = 0.287 ±0.020 and ^^(lO GeV^) = 0.249 ± 
0.015 corresponding to a^iMl^) = 0.118±0.003 [|lO|. As 
Qq = 10 GeV^ is close to the charm threshold, a small 
uncertainty has been included to account for the differ- 
ence between the perturbative QCD corrections for three 
and four flavors. This uncertainty is also included in the 
error estimate for the Bjorken sum rule prediction pre- 
sented in the next section. 
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FIG. 19. Measurements of g\, gf and gi evolved to 
Q'o =5 GeV^ The SMC and E143 g1 data are obtained from 
gl and gf. 



the EUis-Jaffe sum rule prediction. From the combined 
analysis of F^' the EUis-JafFe sum rule is violated by more 
than two standard deviations. The implications of this 
result on t he sp in content of the proton will be discussed 
in Section 



VIII 



VI. RESULTS FOR g^ AND ITS FIRST MOMENT 

A. Evaluation of 32(^1 Q^) 

The spin-dependent structure function ^2 is evaluated 
from the A2 data (Table VIII) using 



92 



2Mx 



7(7 - ri) 
1+7' 



D 



2Mx 



(6.1) 



from Eqs. (2.7) and ( p.9| ) and a parametrization of A\\/D 
from Refs. [||J|. We assume that ^/Q^A^ and ApD 
are independent of Q^ which is consistent with the data. 
The new analyses of g^ or F2 do not affect the (^2 results 
that we published in Ref. |l^ due to the limited accu- 



racy of the data. The gj values are given in Table XIV . 
The expected values of g2^^ and the upper bound of 
§2, based on the positivity limit of A2 are also included. 
The statistical accuracy on 52 is poor since the error is 
proportional to l/x^ and \/Q^, and the data are charac- 
terized by small- a; and large Q^ . All values are consistent 
with zero. 



We re-evaluated the first moments for all experiments 
at their average Q^ using the gi evolution described in 
Section VB. In Figure M the results are shown as a 
function of Q^ . All experimental results are smaller than 
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FIG. 20. Comparison of the experimental results for r!j' to 
the prediction of the EUis-Jafle sum rule. 



B. The first moment of ^j 

The Burkhardt-Cottingham sum r ule p redicts that the 
first moment of g^ vanishes (Section II G ) . This integral 
is evaluated over the measured x range at the mean Q^ 
of the data (Qq = 5 GeV^) by assuming a constant value 
of \/(^A2{x) within each x bin. We obtain 



1.0 < 



gl{x,Ql)Ax<2.l, 



(6.2) 



0.006 



at 90% confidence level. Our measurement of 32 is 
not accurate enough to perform a meaningful extrap- 
olation to X = Q using the expected (72 Regge behav- 
ior, g2{x — > 0)^a;~^+" Q and to test the sum rule. 
The first moment F2(Q§) can be divided into F2(Qo) = 
r2(Qn)^^ + r^(Qn), where F^w is obtained from gfw 
(Eq. (2.4S)) and F2 from the 32 component. Using a 
parametrization of all gX/Ff data |§,§,|| we find that 
the twist-2 part is, as expected, compatible with zero 
(r2((9o)^^ - 0.001 ±0.008.) A violation of the sum rule 
caused by the '§2 term cannot be excluded by the present 
data. 
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TABLE XIV. Results for the spin-dependent structure function gj- The predicted twist-2 term for §2^ 
the upper limit obtained from \A2\ < \fR arc also given. Only statistical errors are shown. 



Eq. (|2.4q)] and 



X range 


{x) 


m (GeV^) 


{y) 


92 


92 


upper 
92 


0.006-0.015 


0.010 


1.36 


0.72 


0.79 ± 75.84 


0.716 ±0.221 


429 ±61 


0.015-0.050 


0.026 


2.66 


0.57 


7.14 ±13.92 


0.447 ±0.069 


101 ± 12 


0.050-0.100 


0.069 


5.27 


0.42 


1.06 ± 4.77 


0.187 ±0.019 


nA±A.& 


0.100-0.150 


0.121 


7.65 


0.34 


-0.95 ± 2.92 


0.037 ±0.015 


6.1 ±2.8 


0.150-0.300 


0.199 


10.86 


0.30 


0.20 ± 1.66 


-0.073 ± 0.007 


1.9 ±1.2 


0.20 -0.600 


0.378 


17.07 


0.25 


0.65 ± 0.64 


-0.096 ± 0.005 


0.2 ±0.5 



VII. EVALUATION OF Tf TJ AND TEST OF 
THE BJORKEN SUM RULE 

We first test the Bjorken sum rule at Ql ~ 10 GeV^ 
assuming 



al 



9i 



9i 



1 



I^D 



(7.1) 



For this test we employ our present proton data and 
our previously published deuteron data lO-O]. For the 
probability of the deuteron to be in a D-state we have 
taken ujd — 0.05 ± 0.01 which covers most of the pub- 
lishe d va lues [106|. Using the method described in Sec- 
tion VD to account for the correlations between errors 



we obtain 



n 



ri = 0.183 ± 0.034 (gg = 10 GeV2 ) , (7.2) 



where statistical and systematic errors are combined in 
quadrature. The theoretical prediction at the same Q^, 
including perturbative QCD correctio ns up to 0{a^) and 
assuming three quark flavors (Section HE 1 ), is F^* — T" = 
0.186 ±0.002. 

We have also performed a combined analysis of all pro- 
ton and deuteron data at Qg = 5GeV^ (Fig. |l9|). The 
combined Ff is obtained using the same method as de- 
scribed in Section VD for F^'. We find 



FP-F5' = 0.199 ±0.025 (Qg^SGeV^, 
All proton and deuteron data) . 



(7.3) 



r? 



The corresponding theoretical expectation is F^' 

0.181 ± 0.003, which agrees with the experimental result 

as shown in Fig. |2l| . 

The structure function gf of the neutron has also been 
measured by scattering polarized electrons on a polar- 
ized ^He ta rget J s]. The re-analyzed neutron data on gf 
from E142 [107| are included in the combined analysis. 
This requires the combination of F^ , F? and Tf via a fit 
constrained by the integral of Eq. (7.1) and the use of 
a Monte Carlo method to compute the 3x3 correlation 
matrix between F^', F" and Ff . The F^' and Ff are ob- 
tained as before; F" is obtained from the E142 data in 
their measured region, but the small- a; extrapolation is 
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FIG. 21. Comparison of the combined experimental results 
for Fj, r" and F^ with the predictions for the Bjorken and 
the Ellis-Jaffe sum rules. The Ellis-Jaffe prediction is shown 
by the black ellipse inside the Bjorken sum rule band. 



determined from the g" values obtained from the SMC 
proton and deuteron data. The result is 



r? 



F5' = 0.202 ±0.022 (Qg^S GeV^, 
All proton, deuteron and neutron data) . (7.4) 



As discussed in Ref. ]108( | , the central value and the error 
of F" is very sensitive to the SMC proton and deuteron 
data. 

The relation between F?, Ff , F" and the Bjorken sum 
rule is illus trated in Fig. |2l| and the results are given in 
Tabic XIIF Proton, deuteron and neutron results con- 
firm the Bjorken sum but disagree with the Ellis-Jaffe 
sum rule. 
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VIII. SPIN STRUCTURE OF THE PROTON 



A. The X dependence of g^ and g^ 

In Figure O we show our results for g\ and g^, together 
with g^ obtained from g\ and gf using Eq. ([7\l|). We 
find that the ratio gi/g\ is close to -1 at small-x, in 
contrast to the r atio F2/F2 which is close to +1 for a; < 
0.01 [pl p09|pl^ . In the QPM the difference between gf 
and gf can be written as 



9l 



5" = -[Awv(a;) - Adv(x) + 2Au{x) - 2Ad{x)] 



(8.1) 



Under the assumption of flavor symmetry in the po- 
larized sea (Au ~ Ad) 111 112 ], the small- a; behav- 
ior of gf/gf indicates a dominant contribution from 
the valence quarks. This is consistent with our results 



from semi-inclusive spin asymmetries |14 which show 
that [Auv{x) — Adv{x)] is positive and that Au^^x) and 
Adv{x) have opposite signs. Fits of polarized parton dis- 
tribu tions in the NLO analysis lead to the same conclu- 
sion poc|,pm|]. 



B. The axial quark charges 

When only three flavors contribute to the nucleon spin, 
the first moment of g^ can be expressed in terms of the 
prot on ma trix elements of the axial vector currents (Sec- 
tion |[I E l| ) 



r?(Q^) = 



cfHQ') 



12 



1 

0-3 + o"8 



C!iQ' 



-ao(Q') 



(8.2) 



We obtain ao{Q^) from r^{Q^) and the experimental 
non-singlet matrix elements 03 and og, which are calcu- 
lated from 5a/<?v and F/D, as presented in Section |V B| . 
The singlet (non-singlet) coefficient function Cf (Cf ) is 
the same as presented in Section [I E 1 , and Cf is com- 



puted with the coefficients c^ in the last column of Ta- 
ble 0. If instead the coefficient from the third column 
were used, we would get ag°. Numerically, a"^ is smaller 
than ao{Q^) by 10% at Q^ = lOGeV^. 

From the combined analysis of all proton data we find 



ao(Q^):= 0.37 ±0.11 



(Q2 ^ 5 GeV2, 
All proton data ) . 



(8.3) 



In Table XV we compare the results with those b ased 
on SMC data only. Calculations in lattice QCD |ll3| 
agree with the measured values of both oq and g^l 9v- 
Using 03 = Ou — fld, as = flu + Od — 2as and ao(Q^) = 
ttu + Id + Og, the individual contributions from quark 
flavors are evaluated from 



Ou = ^ [2ao(Q^) + as -f Sas] , 
Od = - [2ao((5^) + as - Sas] , 
as = - [ao(Q^)-as] . 



(8.4) 
(8.5) 
(8.6) 



The results are given in Table XV. They indicate that as 
is negative, in a greemen t with the measurement of elastic 
v-Y> scattering |114, 115| . 

In the QPM, Oj — Aq^. However, as discussed in Sec- 
tion [IF, due to the U(l) anomaly of the singlet axial 



vector current the axial charges receive a gluon contri- 
bution. In the AB scheme p^ used in our QCD fit for 
three fiavors we have 



A?. - ^4?^A5(Q') {^^u,d,s) 



(8.7) 



In this scheme Aqi is independent of Q^ . For this reason 
some authors consider this to be the correct scheme when 
assuming As = [|3-|9). 

The relation between the matrix element a^ and the 
neutron /3-decay constant gA/gv relies only on the as- 
sumption of isospin invariance. However, in order to re- 
late as to the semileptonic hyperon decay constants F 
and D, we assume SU(3) flavor symmetry and hence con- 
clusions on ao depend on its validity. SU(3) symmetry 
breaking effects do n ot va nish at first order for axial vec- 
tor matrix elements | 116 ], as they do for vector matrix 
elements [I17|. It has been suggested that in order to 
reproduce the experimental values of F and D, the QPM 
requires large relativistic corrections which depend on the 
quark masses; since the s quark mass is much larger than 
that of u and d quarks, these corrections should break 
SU(3) symmetry. Similarly, the relations between the 
baryon mag netic moments predicted by SU(3) are badly 
broken [TT^ . 

The uncertainty on as propagates into gq and as ac- 
cording to 



dap 
das 
dos 
das 



-0.23 



12Cf 



-0.44 



(8.8) 



(8.9) 



The smaller magnitude of Ug and its larger derivative with 
respect to as make it much more sensitive to uncertainties 
in as than ao [HE]. For instance, the experimental test 
of SU(3) from the compatibility of different hyperon /3 
decays allows for a 15% modification of as; this would 
change Og by as much as 50%, while ao changes by less 
than 10%. 

A result for og has been obtained from a leading-order 
l/Nc expansion ]120| ] which is much smaller than the 
value based on the SU(3) analysis. The use of this smaller 
value of og causes og to become larger, while Os becomes 
positive. 
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TABLE XV. Results for ao and individual quark contributions from proton data. The results based on SMC data only are 
given at the average Q^ of the data, Qq — 10 GeV^, and at Qq = 5 GeV"^ for a direct comparison with the combined analysis 
of all proton data. 



Data used 


ao 


flu 


ad 


fls 


SMC rP(10 GeV^) 


0.28 ±0.16 


0.82 ± 0.05 


-0.44 ±0.05 


-0.10 ±0.05 


SMC rP(5 GeV^) 


0.28 ±0.17 


0.82 ± 0.06 


-0.44 ± 0.06 


-0.10 ±0.06 


All rP(5 GeV^) 


0.37 ±0.11 


0.85 ± 0.04 


-0.41 ±0.04 


-0.07 ±0.04 



In principle another source of uncertainty arises from 
the possible contributions of heavier quarks. The heavy 
quark axial current has a non-zero matrix element be- 
cause it can mix with light quark operators |121]. This 
mixing is closely related to the U( l) anomaly and is 
directly calculable in QCD |122,123|, where the heavy 
quark contributions can be expressed in terms of light 
quark contributions. Following th e an alysis of Ref. |121| 
and using the result for ao of Eq. (3^), the expected val- 
ues for flb and a^ for Q^ ^ m\ are —0.003 ± 0.001 and 
—0.006 ± 0.002, respectively. In view of the current ac- 
curacy for ao and of the Q^ range covered by the data, 
the contribution from heavier quarks can be neglected. 



C. The spin content of the proton 



The nucleon spin can be written: 



2 



Lq ± A.g ± Lg 



(8.10) 



Lg) are 



in which AS = Au ± Ad + As and Ag are the con- 
tributions of the quark and gluon spins to the nucleon 
spin, and Lq and Lg are the components of the orbital 
angular momentum of the quarks and the gluons along 
the quantization axis |124]. The Q^ dependence of the 
angular momentum terms analyzed in LO was studied 
in Ref. [125|. It is observed that the asymptotic limit 
(Q^ -^ oo) of the terms (^AS + Lg) and (Ag 
about the same and equal to ^ 1/4. 

In the naive QPM, Ag = Lg — Q and AE = ao. In 
this framework earlier experiments concluded that only 
a small fraction of the nucleon spin is carried by the quark 
spins and that the strange quark spin contribution is neg- 
ative. This conclusion is in disagreement with the Ellis- 
Jaffe assumption of As = Og = 0, which corresponds to 
AE = ag ~ 0.57 with Lq carrying about half of the to- 
tal angular momentum. The Skyrme Model also assumes 
Ag — Lg = 0. In a recent version of this model, where 
9a/ gY is calculated to within 4% of the experimental 
value, AE is found to be between 0.18 and 0.32 [126|. 

In QCD ao differs from AE in a scheme dependent 
way. In the AB scheme the determination of AE and the 
various Ag^ from the measured ao and ai requires an in- 
put value for /S.g. The allowed values for AE and f or t he 
Agi are shown in Fig. p3 as a function of Ag (Eq. (^.7D). 




Ag 



FIG. 22. Quark spin contributions to the proton spin as 
a function of the gluon contribution at Q^ = 5 GeV^ in the 
Adler-Bardeen scheme. 



We see that a value of As — and AE ~ 0.57 corre- 
sponds to Ag{Q^) « 2 at Qo = 5GeV^. However, the 
gluon contribution Ag could be smaller than indicated in 
Eq. ( p77 ) due to finite quark masses and a possibly non- 
negligible contribution from charm, according to the au- 
thors of Ref. |127|. In the absence of direct measurements 
of Ag our results can only be compared with the estimate 
of Ag{Q^) obtained from NLO GLAP fits to the gi data 
Different estimates of Ag{Q^) have 



as in Section VB. 
been obtained. The facto rization scheme used in the fit of 
Ref. |l| and Section [V^ provides AE and Ag{Q'^), while 
ao((3^and Ag{ Q^) a re obtained in the scheme used for 
the fits of Refs. (lO|] and (lO^. While the singlet dis- 
tribution depends on the factorization scheme, the gluon 
distribution is the same in both |5^. Fo r Qn = 5 GeV^ 
we find Ag{Ql) = 1.7 ± 1.1 and Refs. |28| and |lOC| ] 
find Aff(Q n) — 2.6 and 0.76, respectively. The results of 
Ref. [[l28[ are based on the method of Ref. [Q. Simi- 
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larly, at Qq = lOGeV^ it is found that A5i(Qq) is equal 
to 2.0 ± 1.3, 3.1, and 0.89, respectively. 



D. Combined analysis of ao from all proton, neutron 
and deuteron data 



The analysis used to test the Bjorken sum rule can be 
extended to evaluate hq, giving 
( Proton, deuteron and neutron data, Qq = 5 GeV^) , 



ao= 0.29 ±0.06 
fld = -0.43 ±0.02 



Ou = 0.82 ±0.02, 
a, = -0.10 ±0.02 . 



An analysis of oq based on a different selection and 
treatment of experimental data has been presented in 
Ref. [129 1, with similar results. 



IX. CONCLUSION 
A. Summary 

We have presented a complete analysis of our measure- 
ment of the spin-dependent structure function gi of the 
proton from deep-inelastic scattering of high-energy po- 
larized muons on a polarized target. The data cover the 
kinematic range 0.003 < x < 0.7 for Q'^ > I GeV^, with 
an average Q^ = 10 GeV . In addition to these data we 
have also shown for the first time virtual photon-proton 
asymmetries in the kinematic range 0.0008 < x < 0.003 
and Q^ > 0.2 GeV^. In the kinematic range x < 0.03 
our data are the only available measurements of the spin- 
dependent asymmetries. 

The virtual photon asymmetry A^ ~ 9il ^i shows no 
Q^ dependence over the x range of our data within the 
experimental uncertainty. This observation holds when 
we combine our results with those from electron scatter- 
ing experiments performed at smaller Q^ . However, gi 
and Fi are predicted to evolve differently and the dif- 
ference should be observable at small- a; in future precise 
measurements. 

From our data on g\ together with our deuteron data 
we find that the ratio gi/gf is close to -1 at small-x (~ 
0.005), in contrast to F2/F2 which approaches +1. This 
suggests that either the valence quarks give a significant 
contribution to the net quark polarization in this region, 
or that the spin distribution functions of the u and d 
sea quarks are different, i.e. Au{x) ^ Ad{x). The data 
suggest a rise in gi{x) as x decreases from 0.03 to 0.0008. 
A low- a; extrapolation of gi beyond the measured region 
is necessary to compute its first moment Fi and test sum 
rule predictions. Precise data at low x are crucial for 
constraining this extrapolation. 

The new data have initiated much theoretical activ- 
ity in recent years, resulting in an extensive discussion 
of the NLO QCD analyses of the x and Q^ dependence 



of gi and of the interpretation of ao in terms of the spin 
content of the nucleon. As a result, we have used new 
methods for the evaluation of the structure function gi at 
fixed Q^. From this evolved structure function we deter- 
mined the first moment of gi and confirmed the violation 
of the EUis-Jaffe sum rule for the proton by more than 
2(7. We obtain for the singlet axial charge of the proton 
aoiQl) = 0.28±0.16 at Qg = lOGeV^ From the fit to all 
currently available data, we obtain Ag((3§) = 2.0 ±1.3 
which in the Adler-Bardeen renormalization scheme im- 
plies the value AE ~ 0.5. The new data and theoretical 
developments now afford a first glimpse of the polarized 
gluon distribution and its first moment. 

The Bjorken sum rule is fundamental and must hold 
in perturbative QCD. When corrections up to 0{al) 
are included, it predicts F^ - F^' = 0.186 ± 0.002 at 
Qq = 10 GeV^. Using the first moments of the struc- 
ture functions gi evaluated from our proton and deuteron 
data, we find F^ - F^' = 0.183 ± 0.034 at Ql = 10 GcV^ 
in excellent agreement with the theoretical prediction. 
Combining our data with all available data results in a 
somewhat more precise confirmation of the Bjorken sum 
rule. 



B. Outlook 

New data on the spin-dependent structure functions 
gi and 172 of the nucleon are expected in the next two 
years from the SMC, the E154 and E155 collaborations at 
SLAG, and from the HERMES collaboration at HERA. 
However, further knowledge is needed of the low x behav- 
ior of gi and of the polarized gluon distribution Ag(x) 
due to the limited coverage in x and Q^ of these experi- 
ments. 

Future experiments are planned at various experimen- 
tal facilities, including semi-inclusi ve p olarized proton- 
proton scattering by RHIG SPIN Jl30t at BNL, semi- 
inclusive polarized muon-nucleon scattering by COM- 
PASS [131 1 at GERN, and a similar semi-in clusi ve po- 
larized electron-nucleon experiment at SLAG [132[. Fur- 
thermore, a polarized electron-proton collider experi- 
ment at HERA to study the inclusive and semi-inclusive 
scattering is also under consideration |133{ . The non- 
Regge behavior of the unpolarized structure function F2 
has been observed at H ERA in agreement with perturba- 
tive QCD predictions [134, 135 [. The corresponding be- 
havior predicted for the polarized spin structure func- 
tion gi is particularly interesting due to the fact that 
the higher-order corrections in the polarized case are ex- 
pected to be stronger [12? ,13(:[. Also, unlike the unpolar- 
ized case where only the gluon distribution is important 
at low X, in the polarized case the singlet quark, the non- 
singlet quark, and the gluon distributions all play a role. 

In conclusion, the study of the spin structure of the 
nucleon appears certain to remain active well into the 
next century. 
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